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Although the Poincare' and the geometrization conjectures were recently 
proved by Perelman, the proof relies heavily on properties of the Ricci flow pre- 
viously investigated in great detail by Hamilton. Physical realization of such 
' a flow can be found, for instance, in the work by Friedan (Ann. Physics 163 

. (1985), 318-419). In his work the renormalization group flow for a nonlinear 

sigma model in 2 + e dimensions was obtained and studied. For e = 0, by ap- 
proximating the /3— function for such a flow by the lowest order terms in the 
sigma model coupling constant, the equations for Ricci flow are obtained. In 
^ , view of such an approximation, the existence of this type of flow in nature is 

' questionable. In this work, we find totally independent justification for the exis- 

tence of Ricci flows in nature. This is achieved by developing a new formalism 
"^h ' extending the results of two dimensional conformal field theories (CFT's) to 

, three and higher dimensions. Equations describing critical dynamics of these 

O 1 CFT's are examples of the Yamabe and Ricci flows realizable in nature. Al- 

ii , though in the original works by Perelman some physically motivated arguments 

' can be found, their role in his proof remain rather obscure. In this paper, steps 

, are made toward making these arguments more explicit thus creating an op- 



portunity for developing alternative, more physically motivated, proofs of the 
Poincare' and geometrization conjectures. 
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1 Introduction 



The history of physics is full of examples of situations when experimen- 
tal observations lead to very deep mathematical results. Beginning with New- 
ton's dynamics equations, which gave birth to analysis, through the heat and 
wave equations, which gave birth to mathematical physics, through Maxwell's 
equations (experimentally discovered by Faraday), which ultimately had led to 
special relativity on one hand and the theory of homology and cohomology on 
another, etc. thus stimulating the development of topology. Methods of ce- 
lestial mechanics, especially those developed by Poincare', had led Bohr and 
his collaborators to quantum mechanics, on one hand, and to KAM and chaos 
theory, on another. A simple observation by Einstein that gravity can be lo- 
cally eliminated in the freely falling elevator cabin had led to general relativity, 
on one hand, and to great advancements in differential geometry, on another. 
Based on these examples, it is only natural to expect that the recent complete 
proof of the Poincare' and geometrization conjectures by Perelman [1-3] should 
be connected with some processes taking place in the real world. The purpose 
of this work is to demonstrate that this is indeed the case. Since this paper 
is not a review, we made no efforts to list all the latest physics-related results 
inspired by Perelman's works. Instead, we shall only quote those works which 
are immediately relevant to the content of our paper. 

We begin with the observation that in the work of Friedan [4] on the non- 
linear sigma model in 2 + e dimensions, published in 1985 1 , the reference was 
made to a talk by Bourguignon (delivered in 1979) "On Ricci curvature and 
Einstcn metric" in which for the first time the idea of usefulness of Ricci flows 
to differential geometry and related disciplines was introduced, e.g. see Rcf.10 
in Friedan's paper. This idea was picked up by Hamilton in 1982 [5] who began 
from that time on a systematic study of Ricci flows. Although written with 
exceptional clarity, his results summarized in Ref. [6] apparently escaped physi- 
cists attention for many years. This is caused by several factors, which are worth 
mentioning. The Ricci-type flow equations appeared in the physics literature 
in the 80's in connection with the renormalization group analysis of nonlinear 
sigma models and also of strings propagating on curved backgrounds [4,7] with- 
out any reference to the parallel developments in mathematics. Because of this, 
the point of view was developed in the physics literature, especially that related 
to string theory [8], that all major results of general relativity along with quan- 
tum corrections can be deduced from the action functional whose minimization 
produces the Euler-Lagrange equations for the renormalization group (RG) flow 
for strings in curved backgrounds, e.g. see [8], Ch-r 3, page 180. In the simplest 
case such a functional ( e.g. see Eq.(3.4.58) on page 180) is a modification of the 
Hilbert-Einsten functional describing gravity, leading to the action functional 

x It is based on his PhD thesis completed in August of 1980. 
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for the dilaton gravity. The only problem with such a functional lies in the 
fact that it is "living" in 26 dimensional space-time for a bosonic string and 
in 10 dimensional space-time for a fermionic (heterotic) string [9]. These facts, 
apparently were known to Grisha Perelman who used them in his calculations, 
many of which (but not those involving such functionals) are specific only for 
3 dimensional Riemannian manifolds and, hence, are not readily extendable to 
higher dimensions and a pscudo-Ricmannian-type of metric. 

At this point it is essential to notice that both in the work by Friedan [4] 
and in the string-theoretic literature just cited, the dilaton gravity functional 
emerges only as a result of the lowest order approximation in RG-type calcula- 
tions. This fact can be seen very clearly on page 322 of Fricdan's paper, where 
in his Eq.s (1.5) and (1.6) one has to put e = and retain only terms of or- 
der 0{T ) in the sigma model coupling constant T _1 . Such a retention then 
straightforwardly leads to the equation for the Ricci flow. 2 When terms of or- 
der T and higher are included into these RG flow equations, such an RG flow 
is no longer of Ricci type. The procedure of an ad hoc truncation of an infinite 
asymptotic RG series expansion for the RG (3— function (e.g. see Eq.(1.5) in 
Friedan's paper) in order to get the needed equation for the Ricci flow makes 
the physical existence of such a flow doubtful. Since both the Poincare' and the 
geometrization conjectures were proved with help of the unperturbed (nonmod- 
ified by terms other than those induced by diffeomorphisms) Ricci flows [1-3], 
it is essential to find physical processes that employ such unperturbed Ricci- 
type flows. The present work provides examples of such physical processes. It 
enables us to reconnect recent significant mathematical advances with physi- 
cal processes realizable in nature in accord with the logic of development of 
mathematical physics described at the beginning of this section. 

In view of space limitation, before discussing the content of this paper in 
some detail, we would like to mention the following. The existing proofs of the 
Poincare' and the geometrization conjectures are either too compressed [1] or 
too bulky [2,3,10]. In both cases, it may take a considerable amount of time 
for the nonexpert reader to learn and to understand these papers. Fortunately, 
there are succinct publications the present the key ideas of the proof in a lan- 
guage that is considerably more familiar to a much wider group of people with 
mathematical background. Among these reviews, we would especially recom- 
mend those written by Bessieres [11], on the proof of the Poincare' conjecture, 
and by Kapovich [12] on the proof of the geometrization conjecture. This pa- 
per is written under the assumption that our readers have previously read these 
reviews. 

The organization of the rest of this paper is the following. In Section 2 
the scaling analysis of the Ginzburg-Landau (G-L) functional, widely used in 
the theory of critical phenomena, is discussed using both arguments known in 

2 At the same time, to obtain the generating functional whose minimization produces the 
Euler-Lagrangc cquaton for such a flow is nontrivial. In mathematics literature the credit for 
finding such a functional goes to Perelman. In physics literature this functional, Eq. (3.4.58) 
on page 180 of Ref. [8], was known much earlier. It describes the dilaton gravity (albeit in the 
nonphysical space-time dimensions). 
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physics literature and those known in mathematics. In particular, we argue that 
in the physics literature there is a marked difference in treatments of two and 
the higher dimensional models exhibiting critical behavior. While in two dimen- 
sions the full conformal invariance is used (at least at criticality), in dimensions 
higher than two only scale invariance is taken into account at criticality. Using 
scaling arguments from both physics and mathematics, we demonstrate how the 
existing scaling treatments in higher dimensions can be improved in order to 
take into account the effects of full conformal invariance. This improvement is 
brought to completion in Section 3 where we connect the G-L functional used in 
the physics literature with the Yamabe functional known in mathematics. We 
demonstrate that, if properly interpreted, both functionals upon minimization 
produce the same G-L-type equations used in literature on critical phenomena. 
The advantage of working with Yamabe functionals lies in our observations 
that a) such a functional can be rewritten in terms of the Hilbcrt-Einstein 
action functional for Euclidean gravity and, b) that such a functional is man- 
ifestly conformally invariant. Such a connection between the H-E action and 
its pseudo-Riemannian extension is needed in applications in relativity which 
are studied in Section 4. Connections between the G-L and Yamabe functionals 
exist only for dimensionality higher than two. In two dimensions one has to find 
an analog of the G-L- Yamabe functional. This task is accomplished in Section 
5 in which not only such an analog is found but, also, its connection with known 
results in string and CFT. Specifically, we were able to connect such an analog 
with the much earlier results obtained by Distler, Kawai and David for 2d quan- 
tum gravity in noncritical dimensions pedagogically summarized in the book by 
Hatfield [13]. Applications of their results to CFT are summarized in the lec- 
ture notes by Abdalla et al [14] and, more recently, in the review by Nakayama 
[15]. Because of such connections, the extension of these 2 dimensional CFT 
results to higher dimensions can be done quite systematically. This is discussed 
in Section 6. In Section 7 we develop an alternative method allowing us to 
reobtain results known already for 2 and higher dimensions. This formalism is 
based on the observation that results obtained in previous sections are just static 
solutions of more general equations describing critical dynamics for systems of 
the G-L-type. By noticing that mathematically such dynamics coincides with 
that known for the Yamabe-type flows discovered by Hamilton [16], we follow 
the logic of his works in order to relate the Yamabe flow to the more general 
Ricci flow of which the Yamabe flow is a special case. To do so we use some 
results by Perelman [1], which he was using in his proofs of the Poincare' and 
the geometrization conjectures. Our use of Perelman's results is not mechanical 
however. By exposing physical arguments hidden in his work, we were able to 
reobtain many of his results much simpler as compared to treatments which 
can be found in current works by mathematicians [2,3,10] aimed at explaining 
and elaborating on many hidden details of Perelman's original proof. Finally, in 
Section 8 we discuss other possible physical systems whose dynamics is expected 
to be described in terms of Ricci flows. We also discuss the conditions under 
which experiments on critical dynamics should be conducted in order to detect 
the effects of Ricci and Yamabe flows. 
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2 Scaling analysis of the G-L functional 



Conventional scaling analysis of the G-L functional used routinely in physics 
literature can be found, for example, in Ref.[17]. This analysis differs somewhat 
from that for the </> 4 model as described in the monograph by Itzykson and 
Zubcr [18]. For the sake of the discussion that will follow, we would like to 
provide a sketch of the arguments for both cases now. 

We begin with the <f> model following Ref. [18]. Let C(x) be the Lagrangian 
of this scalar field model whose action functional in d dimensions S[(f>] is given 
by S[4>] = J d d xC{x). Let furthermore A be some nonnegative parameter. Then 
the requirement that S[4>\ be independent of A, i.e. J d d xC{x) = J d d x\ L(Xx), 
leads to the constraint 



/ 



d d x(x ■ 4- + d)C(x) = (2.1) 
ox 



obtained by differentiation of S[4>] with respect to A with A being equalled to 
one at the end of calculation. For C(x) given by 

A*) = ^(V*) 2 + ^ 2 + ^ 4 (2-2) 
the change of C(x) under the infinitesimal scale transformation is given by 3 

g = (x ■ ^ + d)C(x) +(d 4)g</> 4 - to 2 2 . (2.3) 

Comparison between Eq.s(2.1) and (2.3) implies that the action S[4>] is scale 
invariant if d = 4 and m 2 = 0. The result just obtained raises an immediate 
question: given that the massless G-L action is scale invariant for d = 4, will 
it also be conformally invariant under the same conditions? We provide the 
answer to this question in several steps. 

First, let M be some Riemannian manifold whose metric is g. Then, any 
metric g conformal to g can be written as g = e^g with / being a smooth real 
valued function on M [19]. Let A g be the Laplacian associated with metric g A 
and, accordingly, let Ag be the Laplacian associated with metric g. Richardson 
[20] demonstrated that 



A s = e-fA g -^-l)fe-fA g -^-l)e-f(A g f) + ^-l)e-fA g of (2.4) 



(here (A g o /) * should be understood as A g (f(x)^/(x)). Eq.s(2.3) and (2.4) im- 
ply that the conformally invariant (string-type) functional S[X] = J M d 2 x^/g(\/ g X^)- 



3 In arriving at this result we took into account Eq. (13-40) of Rcf. [18] along with condition 
D = | - 1, with D defined in Eq.(13-38) of the same reference. Also, we changed signs (as 
compared to the original source) in C(x) to be in accord with the accepted conventions for 
the G-L functional. 

4 That is Ag* = - (det g) ~ 5 d t (g ij (det g) 5 dj *) for some scalar function *(x). 
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(VgX) exists only for d = 2. For d > 2, in view of Eq.(2.4), the conformal in- 
variance of the </> 4 model is destroyed. Thus, using conventional field-theoretic 
perturbational methods one encounters a problem with conformal invariance at 
the zeroth order level (in the coupling constant G), unless d — 2. This problem 
formally does not occur if one requires our physical model to be only scale 
invariant. This is undesirable however in view of the fact that in 2 dimensions 
an arbitrary conformal transformation of the metric tensor g of the underlying 
two dimensional manifold M is permissible at criticality [21]. Abandoning the 
requirement of general conformal invariance in 2 dimensions in favor of scale 
invariance in higher dimensions would destroy all known string-theoretic meth- 
ods of obtaining exact results in two dimensions. Our general understanding 
of critical phenomena depends crucially on our ability to solve two dimensional 
models exactly. All critical properties for the same type of models in higher 
dimensions are expected to hold even in the absence of exact solvability. Fortu- 
nately, the situation can be considerably improved by reanalyzing and properly 
reinterpreting the scaling results for dimensions higher than 2. 

This observation leads to the next step in our arguments. Using the book 
by Amit [17], we consider first the scaling of the non interacting ( free ) G-L 
theory whose action functional is given by 



S[</>] = J d d x{(v4>) 2 +m 2 2 }. (2.5) 

Suppose now that, upon rcscaling, the field <p transforms according to the rule 5 : 

~(f>{Lx) = L u ^>(x). (2.6) 

If we require 

J d d x{{sj4>? + m 2 <t> 2 } = j d d xL d {(^) 2 + m 2 2 } (2.7) 

and use Eq.(2.6) we obtain, 

S[<j>] = J d d xL d {(^(f>) 2 L 2uJ - 2 + rrS^L 2 "}. (2.8) 

In order for the functional S[<fr] to be scale invariant the mass m 2 should scale 
as follows: rh 2 — m 2 L~ 2 . With this requirement the exponent u) is found to be: 
w = l-i 

The above scaling can be done a bit differently following the same Ref.[17]. 
For this purpose we notice that although the action S[4>] is scale invariant, there 
is some freedom of choice for the dimensionality of the field <f). For instance, 
instead of S[4>] we can consider 

S[4>} = ^Jd d x{( V ^ 2 +m 2 <f> 2 } (2.9) 



5 Notice that this is already a special kind of conformal transformation 



G 



with a d being some volume, say, a d = J d d x . Then, by repeating arguments 
related to Eq.(2.8), we obtain: u = 1 (instead of w = 1 — | previously ob- 
tained ) in accord with Eq.(2-66) of the book by Amit [17]. Although, from 
the point of view of scaling analysis both results are actually equivalent, they 
become quite different if we want to extend such scaling analysis by considering 
general conformal transformations. Even though, in view of Eq.(2.4), such a 
task seems impossible to accomplish, fortunately, this is not true as we would 
like to demonstrate. This leads us to the next step. 

Noticing that the mass term scales as the scalar curvature R for some 
Riemannian manifold M, i.e. the scaling m 2 = m 2 L~ 2 is exactly the same as 
the scaling of R given by 

R = L- 2 R. (2.10) 

This result can be found, for example, in the book by Wald, Ref. [22], Eq.(D.9), 
page 446. In general, the scalar curvature R(g) changes under the conformal 
transformation g — e 2 -' g according to the rule [19] 

R(g) = e- 2 f{R(g) - 2(d - l)A g f - (d - l)(d - 2) | Vfl /| 2 } (2.11) 

where A g f is the Laplacian of / and Vs/ is the covariant derivative defined 
with respect to the metric g. From here we see that, indeed, for constant / 
's the scaling takes place in accord with Eq.(2.10). Now, however, we can do 
more. Following Lee and Parker [19], we simplify the above expression for R. 
For this purpose we introduce a substitution: e 2i — tp p ~ 2 , where p = -j^, so 
that g — <f p ~ 2 g- With such a substitution, Eq.(2.11) acquires the following form: 

R(g) = <p 1 -r{aA g <p + R(g)<p), (2.12) 

with a — 4^5^ • Clearly, such an expression makes sense only for d > 3 and 
breaks down for d = 2. But we know already the action 5[X] which is both 
scale and conformally invariant in d = 2. It is given after Eq.(2.4). Fortunately, 
the results obtained in this section can be systematically extended in order to 
obtain actions, which are both scale and conformally invariant in 3 dimensions 
and higher. This is described in the next section. 



3 G-L functional and the Yamabe problem 

We begin with the following observation. Let R(g) in Eq.(2.11) be some constant 
(that this is indeed the case, we shall demonstrate shortly). Then Eq.(2.12) 
acquires the following form 

aA g <p + R( g )<p = R{g)<fP- 1 . (3.1) 

By noticing that p — 1 = ^| we obtain at once : p — 1 = 3 (for d = 4) 
and p — 1 = 5 ( for d = 3). These are familiar Ginzburg-Landau values for 
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the exponents of interaction terms for critical and tricritical G-L theories [23]. 
Once we recognize these facts, the action functional producing the G-L-type 
Eq.(3.1) can be readily constructed. For this purpose it is sufficient to rewrite 
Eq.(2.9) in a manifestly covariant form. We obtain, 

S[<p] = -r f d d x^-g{a(x7 gV f + R(g) V 2 } ee f^l (3.2) 

Minimization of this functional produces the following Euler-Lagrange equation 

aAgip + R{g)ip ~ XipP- 1 = (3.3) 
with constant A denoting the extremum value for the ratio: 

A = -j^l = m{{S[ip] : g conformal to g}. (3.4) 

II r II p 

In accord with the Landau theory of phase transitions [24] it is expected that 
the conformal factor ip is a smooth nonnegative function on M achieving its 
extremum value (p . Comparison between Eq.s(3.1) and (3.3) implies that A = 
R(g) as required. These results belong to Yamabe, who obtained the Euler- 
Lagrange G-L-type Eq.(3.3) upon minimization of the functional S[ip] without 
any knowledge of the Landau theory. The constant A is known in literature as 
the Yamabe invariant [19,25]. Its value is an invariant of the conformal class 
(M,g). Given these facts, we obtain the following: 

Definition 3.1. The Yamabe problem lies in finding a compact Riemannian 
manifold (M , g) of dimension n > 3 whose metric is conformal to the metric g 
of a constant scalar curvature. 

Subsequent developments, e.g. that given in Ref. [26,27], extended this prob- 
lem to manifolds with boundaries and to non- compact manifolds. It is not too 
difficult to prove that the (Yamabe-Ginzburg-Landau-likc) functional is mani- 
festly conformally invariant. For this purpose, we need to rewrite Eq.(2.12) in 
the following equivalent form 

^R(g) = (apA g p + R{g)^ 2 ). (3.5) 

This can be used in order to rewrite E[tp] as follows : E[ip] = J d d Xy/gip p R(g). 
Next, by noting that / d d x\fg = J d d x s /gyj p , we can rewrite the Yamabe func- 
tional in the Hilbert-Einstein form 

S[ V ] = S^Sm (3.6) 
(/ cfixy/g) p 

where both the numerator and the denominator are invariant with respect to 
conformal changes in the metric, e.g. scale change. 
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Remark 3.2. In order to use these results in statistical mechanics, we have 
to demonstrate that the extremum of the Yamabe functional S[<p] is realized for 
manifolds M whose scalar curvature R(g) in Eq.(3.3) is also constant. This is 
the essence of the Yamabe problem (as defined above) for physical applications. 
Fortunately, the Yamabe problem has been solved positively by several authors 
starting with Yamabe himself (whose proof contained some mistakes however) 
and culminated in the work by Schocn [28]. Details can be found in the review 
paper by Lee and Parker [19] and in the monograph by Aubin [29]. 

Remark 3.3. In view of the relation J d d xy / g = J d d Xy/g<p p , it is clear that 
for the fixed background metric g Eq.(3.3) can be obtained alternatively using 
the following variational functional 

S[<p] - J d d x^g{a{y gV f + R(g)f 2 } -\J d d x^gyP (3.7) 

where the Lagrange multiplier A 6 is responsible for the volume constraint. Such 
a form of the functional S[tp] brings this higher dimensional result in accord with 
that to be developed below for two dimensions (e.g. see Section 5, Eq.(5.24), 
and Section 7 ). 

Apart from the normalizing denominator, Eq.(3.6) represents the Hilbert- 
Einstein action for pure gravity defined for Riemannian d dimensional space. 
The denominator, the volume V taken to power | , serves the purpose of making 
S[(p] manifestly conformally invariant, Ref.[28], page 150. 

4 From G-L to Hilbert-Einstein functional 



In this section we would like to analyze the significance of the cosmological 
constant term in the Hilbert-Einstein action for gravity from the point of view 
of the G-L model. For this purpose, following Dirac [30] let us consider the 
extended Hilbert-Einstein (H-E) action functional S c (g) for pure gravity with 
an extra (cosmological constant) term defined for some (pseudo) Riemannian 
manifold M (of total space-time dimension d) without a boundary: 

S c (g)= [ R^gd d x + C [ d d x^g. (4.1.) 
Jm Jm 



The (cosmological) constant C is determined based on the following chain of 
arguments. First, we introduce the following. 



Definition 4.1. Let Rij be the Ricci curvature tensor. Then the Einstein 
space is defined as a solution of the following vacuum Einstein equation 

Rij = Xgij (4.2) 

with A being constant. 



From this definition it follows that in both the Riemannian and pseudo- 
Riemannian cases 

R = dX (4.3) 
Following Dirac [30], variation of the action S c (g) produces 

Rij-\gijR+\c 9i j=Q. (4.4) 

The combined use of Eq.s(4.3) and (4.4) produces: 

C = A(^) (4.5) 

With these results, using Eq.s(4.3) and (4.5) we can rewrite Eq.(4.4) as 
follows: 

Rij ~ \gijR +^( d - 2 ) R 9ij = 0- (4.6) 

Remark 4.2. These observations allow us to look at possible pseudo- 
Riemannian extension of the results originally obtained by Yamabe, Refs.[19, 
29, 31], for Riemannian manifolds. Although at first sight this might be of sig- 
nificance importance for high energy physics applications, in a companion pub- 
lication [32] we shall argue that dealing with Riemannian manifolds is already 
quite sufficient so that the Yamabe results can be used in high energy physics 
without changes. It should be clear as well that the earlier introduced notion 
of the Einstein space is applicable to both Riemannian and pseudo- Riemannian 
manifolds. 



In view of this remark, we would like to argue that Eq.(4.6) can also be 
obtained by varying the Yamabe functional, Eq.(3.6). Indeed, following Aubin 
[29] and Schocn [28] , let t be some small parameter labeling the family of metrics: 
gij (*) = 9ij + th-ij . Then, these authors demonstrate that 



dR t 
dt 



= V l V 3 h l] -V 3 V 3 K i -R^h l] (4.7) 



3 '"i ±L '"V 

t=0 



and 

1 




^VMS^,, (4-8) 
t=o z 

where, as usual, |<?| = |det gtj\ . Consider now the Yamabe functional, Eq.(3.6), 

but, this time, written for the family of metrics which belong to the same 
conformal class. We have, 

TZ{g{t)) = {V(t))T f R(g(t))DV(t), (4.9) 
Jm 
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where the volume is given by V(t) = J M d d xy // g(t) and, accordingly, DV(t)~d d x^ 
Using Eq.s(4.7) and (4.8) in Eq.(4.9) and taking into account that the combina- 
tion v l V"* h-ij — V"* Vj K is the total divergence produces the following result: 

d -K(g(t))) = V(0)- 2 -^[[ {Rg^/2-W)h l0 DV(0) f DV(0) 

/ t=0 J M JM 



dt 



~{\ -k I DV(0) f h^DViO)}. (4.10) 

z a JM JM 

If the metric g is the critical point of lZ(g(t)), then 

/ DV{Q) + {\-h{( RDV(0))g l3 =0. (4.11) 

JM L a JM 



From here, multiplication of both sides by g % i and subsequent summation pro- 
duces at once 

R-?-d+{ l --- d )<R>d = Q, (4.12) 

where < R >= J RDV(0) is the average scalar curvature. Eq.(4.12) can 
be rewritten as R =< R > . But this condition is formally equivalent to the 
Einstein condition, Eq.(4.2), in view of Eq.s(4.3), (4.6)! Hence, under such 
circumstances, Eq.s (4.11) and (4.12) are equivalent. 

Remark 4.3. It is important to mention that in general the condition R =< 
R > is more restrictive than the condition, Eq.(4.2), defining Einstein spaces. 
That is to say, the constant Ricci tensor Rij always leads to the constant scalar 
curvature R while the Ricmannian spaces of constant scalar curvature R are not 
necessarily spaces for which the Einstein condition, Eq.(4.2), holds. Hence, they 
may or may not be of Einstein type. This was demonstrated by Herglotz in 1916 
(immediately after Einstein's general relativity was completed) [33, page 148]. 
His results were subsequently generalized in Refs. [34,35] . We shall return to this 
topic in Section 7 where we discuss this problem dynamically following the 
works by Hamilton and Perelman. Using dynamical arguments, we demonstrate 
that the equation R =< R > is compatible with the Einstein condition and, 
hence describes the Einstein spaces. 

In view of this remark, in the rest of this work we shall develop our formalism 
under the assumption that for the set of metrics of fixed conformal class, the 
variational problem for the G-L functional, Eq.(3.7), is equivalent to the varia- 
tional problem for the H-E functional, Eq.(4.1), for pure gravity in the presence 
of the cosmological constant. Because of this noted formal equivalence, presence 
of the cosmological term in the H-E action amounts to volume conservation for 
the G-L variational problem. From here, it should be clear that Eq.(2.12) is 
equivalent to the condition: R =< R > in view of Eq.(3.4). 

The result shown in Eq.(4.12) becomes trivial for d = 2. Physically, however, 
the case d = 2 is important since it is relevant to all known exactly solvable 
models of statistical mechanics treatable by methods of conformal field theories. 
Hence, now we would like to discuss needed modifications of the obtained results 
in order to obtain a two dimensional analog of the G-L theory. 
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5 Ginzburg-Landau-like theory in two dimen- 
sions 



5.1 Designing the two dimensional G-L-Yamabe functional 

From the field-theoretic treatments of the G-L model [36] we know that a 
straightforward analysis based on asymptotic e— expansions from the critical 
dimension (4) to the target dimension (2) is impractical. At the same time, the 
results of CFT and exactly solvable models make sense thus far only in d = 2. 
The question arises: is there an analog of G-L Eq.(3.1) in two dimensions? 
And, if such an analog does exist, what use can be made of such a result? In 
this section we provide affirmative answers to these questions. We demonstrate 
that: a) indeed, such a two dimensional analog of the G-L equation does exist 
and is given by the Liouville Eq.(5.19), b) the functional, Eq.(5.14), whose 
minimization produces such an equation is the exact two dimensional analog 
of the G-L-Yamabe functional, Eq.(3.2), c) these results can be (re)obtaincd 
from the existing string-theoretic formulations of the CFT which were developed 
entirely independently. 

To discuss topics related to items a) and b) just mentioned, we begin with 
the observation that in two dimensions, Eq.(2.4) acquires a very simple form 

A g = e^A g , (5.1) 

where we use a factor of 2 to be in accord with Eq.(2.11) for the scalar curvature. 
According to Eq.(2.11), the scalar curvature in two dimensions transforms like 

R{g) = e- 2 f{R(g) - A g 2f} (5.2) 

while the area dA = d 2 x^/g transforms like 

dA = e~ 2f dA. (5.3) 

These facts immediately suggest that the previously introduced action func- 
tional 

5[X] = / d 2 xV?(v s X) • ( Vs X) (5.4) 

J M 

is conformally invariant. Using results by Polyakov [37] and noting that (VgX) • 
(VgX) = g a/3 d a X^dpX^, we need to consider the following path integral 7 

eM-H9)) = j ' D[4>}eM~l J^x^gg^d^dpcj)). (5.5) 

Here the symbol !F(g) stands for the "free energy" usually defined this way in 
statistical mechanics 8 . Fortunately, this integral was calculated by Polyakov [37] 

7 Without loss of generality, we would like to consider the case of one component field <j> 
only. 

8 Usually, instead of F(g) one writes F{g)/k B T,where T is the temperature and kg is the 
Bolzmann's constant. In the present case the problems we are studying do not require specific 
values for these constants. For this reason they will be suppressed. 
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for two dimensional manifolds M without boundaries and by Alvarez [38] for 
manifolds with boundaries. In this work, to avoid unnecessary complications, 
we shall be concerned with manifolds without boundaries. Although in the 
original work by Polyakov, one can find the final result of calculation of the 
above path integral, the details of this calculation can only be found elsewhere. 
In particular, we shall follow pedagogically written papers by Weisberger [39,40] 
and Osgood, Phillips and Sarnak [41-43] (OPS). 

To begin, let g a p be some reference metric, while g a p is a metric conformally 
related to it, i.e. g a p — exp(— 2if )g a/ 3 . Should the above path integral be 
for the flat (i.e.<? a /3=<5 a/ 3 ) two dimensional manifold, one would have at once 
the result: T=\ lndct' Ao, where the prime indicates that the zero mode is 
omitted. Because the flatness assumption in general is physically incorrect, it is 
appropriate to pose a problem: how is the path integral for the metric g related 
to the path integral for the metric g ? The paper [41] by OPS provides an 
answer, e.g. see Eq.(1.13) of this reference. To connect this equation with that 
known in the physics literature we replace it by the equivalent expression 

]n ^ _ k ^ __^iy iMWM , + R(gm m 

M 

useful in applications to strings and CFT, e.g see Ref.[13], page 637 9 . 

It is worthwhile to provide a few computational details leading to Eq.(5.6). 
Let, therefore, ip i be eigenfunctions of the Laplacian A g with eigenvalues Xi 
arranged in such a way that = Ao < Ai < A2 < • • •, i.e. 

-A g ^ i + A^ = 0. (5.7) 

Then, we construct the zeta function 

00 

c 9 w=Er (5.8) 

1=1 

in such a way that 

dot 'A g = exp(-Cg(0)) (5.9) 
with Cg(0) = ( j- s ( g (s)) s=0 • Using Eq.(5.1), we obtain as well 

—e~ 2v Agipi + Xip i = 0. (5.10) 

In particular, for a constant ip — <p we obtain 

00 



9 In arriving at this result we took into account that the Gaussian curvature K in Ref. [41] 
is related to Ricmannian curvature R as R = 2K . 
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Use of Eq.(5.9) in Eq.(5.11) produces: 



C § (0) = C>) + 2^(^-1). ,r,.;i.2) 



This result was obtained with help of the known fact, Rcf.[41], Eq.(1.9), that 

C 9 (0) = ^ - 1 (5-13) 

with x(M) being Euler's characteristic of a two dimensional manifold M without 
boundaries. In view of the definition, Eq.(5.9), and the Gauus-Bonnet theorem, 
we observe that Eq.(5.6) is reduced to Eq.(5.12) for the case of constant confor- 
mal factor p = ip, as required. 

Using Eq.(5.6) and, following OPS [41], we would like to consider the related 
functional T(p) defined by 

T{p) = \J dA g {\ Vg p\ 2 + - ttx(M) In J dA g e 2 *. (5.14) 

M M 

This functional is the exact analog of the Yamabe functional, Eq.(3.2), in di- 
mensions 3 and higher as we shall demonstrate below and in Section 7. In the 
meantime, in view of Eq.(5.6), it can be rewritten in the following equivalent 
form 10 

T{p) = -67rlndct 'A g +tt(6- x(M))lnA (5.15) 
Let a be some constant, then 

T{p + a)=T{p). (5.16) 

Such invariance signifies the fact that the above action is scale invariant. This 
property is in complete accord with Eq.s(2.7) and (2.8). If we impose a con- 
straint (fix the gauge) : A = 1, then we would end up with the Liouvillc-likc 
action used in CFT. We would like to explain all this in some detail now. 

Following OPS, it is convenient to replace the constraint A = 1 by the 
alternative constraint on the field p: 

pdA g = 0. (5.17) 



M 



To demonstrate that such a constraint is equivalent to the requirement A = 1, 
we note that, provided that the field tp minimizes T(ip) subject to the constraint 
Eq.(5.17), the field 

p = tp - i In J exp(2?P)dA g (5.18) 



M 



10 Here one should understand the word "equivalence" in the sence that both functionals 
produce the same critical metrics upon minimization. 
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minimizes !F(ip) subject to the constraint A = 1. Using these facts, minimization 
of produces the following Liouville equation 

M 

Comparing this result with Eq.(5.2) and taking into account that 

2^ X (M) 



/ exp(2?/>)cL4 9 

M 



R{g) = const (5.20) 



we conclude that, provided that the background metric g is given so that the 
scalar curvature R{g) (not necessarily constant) can be calculated, the Liouville 
Eq.(5.19) is exactly analogous to that previously obtained in Eq.(3.3) (which 
is the same as Eq.(2.12)). In view of this analogy, use of Eq.s.(3.4) and (3.6) 
as well as Eq.s(5.14),(5.15) causes the functional to attain its extremum 

for metric g of constant scalar curvature R(g). To decide if the extremum is 
minimum or maximum we have to consider separately the cases x(M) > 
and x(M) < 0. In the case of x(M) > we have only to consider manifolds 
homcomorphic to S 2 so that x(M) = 2. Fortunately, this case was considered 
in detail by Onofri [44] 11 . Using his work, the following inequality 

In J dA § cxp(V) < J dA g i> + J y dA g IvV-l 2 (5.21) 
s 2 s 2 s 2 

attributed to Aubin [29] and inspired by earlier results by Trudinger and Moser 
[45], is helpful for deciding whether the obtained extremum is a minimum or 
maximum. Here g is the metric of the unit sphere S 2 with a constant Gaussian 
curvature of 1. The metric g conformal to g is given by g = cxp(2-0)g , with tp 
obeying the Liouville equation (just like Eq.(5.19)), where both R(g) and R(g) 
are constant by virtue of the initial choice of g. By combining Eq.s(5.6),(5.14) 
and (5.15) with the inequality (5.21) and taking into account that by design 
lm4g = we obtain, 

~ 37rln f^ = I / ^{IV§^| 2 + 2^}-ln J dA g eM^)>0, (5.22) 

s 2 s 2 

with equality occurring only at the extremum ip = V 1 *, with function tp* being 
a solution of the Liouville Eq.(5.19). It can be shown [44] that: a) such a 
solution involves only Mobius transformations of the sphere S 2 and that, b) the 
functional T(tp) is invariant with respect to such transformations. The case of 
x(M) < is treated in Section 2.2 of the OPS paper, Ref.[41] and leads to the 
same conclusions about the extremality of the functional T{(f) . 



11 Generalization of these two dimensional results to higher dimensional manifolds of even 
dimensions can be found in the paper by Beckner [45] . 
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Corrollary 5.1. The two dimensional results just obtained by design 
are in accord with results obtained in higher dimensions ( discussed in Sections 
3,6 and to be discussed in Section 7). In particular, the functional T{(p) is 
the exact analog of the Yamabe functional S[ip], Eq.(3.2). Since in both cases 
the functionals are "translationally" (actually, "scale") invariant, e.g. compare 
Eq.(3.6) with Eq.(5.16), in both cases, the extremum is realized for the metric 
conformal to the metric of constant scalar curvature, e.g. compare Eq.(3.3) with 
the Liouvillc Eq.(5.19). 

5.2 Connections with string and CFT 



The obtained results allow us now to discuss topic c) listed at the beginning 
of this section. Using known results for string and conformal field theories, 
the results just obtained can be given a statistical mechanical interpretation. 
Following [14,15,46], it is of interest to consider averages of the vertex operators. 
They are defined by 

In \ r n 

( J]exp(/3^(^))) = D[4>] C ^{-S L m\{cMP l ^)), (5.23) 
\i=i / J i=i 

where the Liouville action Sl(<I>) is given (in notations adopted from these ref- 
erences) by 

W) = F-/ d M\v4>\ 2 - Qi?(ff)^ + 87r/2exp(a+^)]. (5.24) 

The actual values and the meaning of the constants Q, jl and a + are explained 
in these references and are of no immediate use for us. Clearly, upon proper 
rescaling, we can bring Sx(0) to the form which agrees with ^((p), defined by 
Eq.(5.14), especially in the trivial case when both x{M) and Ji are zero. When 
they are not zero, the situation in the present case becomes totally analogous 
to that discussed earlier for the Yamabe functonal. In particular, in Section 3 
we noted that the G-L Euler-Lagrange Eq.(3.3) can be obtained either by min- 
imization of the Yamabe functional, Eq.(3.2), (or (4.9)) or by minimization of 
the G-L functional, Eq.(3.7), where the coupling constant A plays the role of the 
Lagrange multiplier keeping track of the volume constraint. In the present case, 
variation of the Liouville action Sl(4>) will produce the Liouville equation, e.g. 
see Eq.s(5.19)-(5.20), which is the two dimensional analog of the G-L Eq.(3.3). 
Such variation is premature, however, since we can reobtain !F(<*p) exactly using 
the path integral, Eq.(5.23). This procedure then will lead us directly to the 
Liouville Eq.(5.19). 

For this purpose, we need to consider the path integral, Eq.(5.23), in the 
absence of sources, i.e. when all (i i = 0. Following ideas of Ref.s [14,15,46], we 
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take into account that: a) 



1 

47T 



dAgR(g) = X (M) = 2-2.g (5.26) 



M 



with g being genus of M and, b) the field <fi can be decomposed into <j> = <j> +<p> in 
such a way that 4> is coordinate-independent and ip is subject to the constraint 
given by Eq.(5.17). Then, use of the identity 

oo 

dxexp(ax) exp(— bexp(jx)) — —b~~T( — ) (5.27) 

7 7 

— oo 

(with r(x) being Euler's gamma function) in the path integral, Eq.(5.23), re- 
quires us to evaluate the following integral 

Q - f 

d(f> exp((f> — x(M))exp((fj, / dA g exp(a + </>)) exp(a + o ) 



M 



r(- s ) 



M 



(At J dAg exp(a + ip)) s (5.28) 



with s being given by 

s = -^-X(M). (5.29) 

Using this result in Eq.(5.23), we obtain the following path integral (up to a 
constant) 

Z[p]= J L»b]exp(-^)) (5.30) 
with functional !F(ip) given by 

fbp) = S L {^- Jx = 0) - 2^"X(M) ln[/i | dA § exp(a +V )]. (5.31) 

M 

This functional (up to rescaling of the field ip) is just the same as the Yamabe- 
like functional T(<p) given by Eq.(5.14). 

Define now the free energy T in the usual way via T = — In Z [ip] (as was done 
after Eq.(5.5)) and consider the saddle point approximation to the functional 
integral Z[(p}. Then, for spherical topology in view of Eq.(5.22), we (re)obtain 
T > with equality obtained exactly when ip — ijj* . Inclusion of the sources 
(or the vertex operators) can be taken into account also, especially in view of 
the results of [43]. 

To understand better the physical significance of the obtained results, it 
is useful to rcobtain them using a somewhat different method. The results 
obtained with this alternative method are also helpful when we shall discuss 
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their higher dimensional analogs in the next section. For this purpose we write 
(up to a normalization constant) 



J cxp{-S L (^)} = J dAe-^ A Z L (A) (5.32) 



where 



Zl(A) = 



J Jm 



dAg exp(a + (f)) - A) exp(S L ((f>; \x = 0)). 



(5.33) 



If, as before, we assume that 



ip, then an elementary integration over 



b produces the following explicit result for Z L {A) : 



Z L {A) = —A" I D\ 
a+ 



M 



dA g enp(a + ip) 



where the exponent u> is given by 



iv = 



X(M)Q 
2a+ 



1. 



exp(-S L (<p;n = 0)) 

(5.34) 

(5.35) 



Finally using Eq.(5.34) in Eq.(5.32) produces back Eq.s(5.30) and (5.31) (again, 
up to a constant factor). An overall "-" sign can be removed by proper normal- 
ization of the path integral. These results can be used for computation of the 
correlation functions of conformal field theories (CFT). Details can be found in 
[14,15,46]. 



6 Designing higher dimensional CFT(s) 
6.1 General remarks 

In previous section we explained a delicate interrelationship between the path 
integrals Eq.(5.5) and (5.23). From the literature on CFT cited earlier it is 
known that, actually, both are being used for the design of different CFT mod- 
els in 2 dimensions. For instance, if one entirely ignores the effects of curvature 
in Eq.(5.5), then one ends up with the Gaussian -type path integral whose cal- 
culation for a flat torus is discussed in detail in Ref.[21], pages 340-343, and, 
by different methods, in our Ref. [47]. By making appropriate changes of the 
boundary conditions in such a path integral (or, equivalently, by considering 
appropriately chosen linear combinations of modular invariants) it is possible to 
build partition functions for all existing CFT models. For the same purpose one 
can use the path integral given by Eq.(5.23) but the calculation proceeds differ- 
ently as explained in Section 5. Since in 2 dimensions the conformal invariance 
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is crucial in obtaining exact results, use of the Gaussian- type path integrals is, 
strictly speaking, not permissible. Fortunately, the saddle point -type calcu- 
lations made for the path integral, Eq.(5.23), produce the same results since 
the extremal metrics happen to be flat [47]. If one does not neglect curvature 
effects in Eq.(5.5), one eventually ends up with the integrand of the path inte- 
gral, Eq.(5.23). This result is a consequence of Eq.(5.6) known as the conformal 
anomaly. If one would like to proceed in an analogous fashion in dimensions 
higher than two one should be aware that there is a profound difference be- 
tween calculations done in odd and even dimensions. We would like to explain 
this circumstance in some detail now. 

In 2 dimensions the conformal invariance of the action, Eq.(5.4), has been 
assured by the transformational properties of the 2 dimensional Laplacian given 
by Eq.(5.1). In higher dimensions, the Laplacian is transformed according to 
Eq.(2.4), so that even the simplest Gaussian model is not conformally invariant! 
This observation makes use of traditional string -theoretic methods in higher 
dimensions problematic. In two dimensions these are based on a two stage 
process: first one calculates the path integral, Eq.(5.5), exactly and, second, 
one uses the result of such a calculation (the conformal anomaly) as an input 
in another path integral, e.g. Eq.(5.23), which is obtained by integrating this 
input over all members of the conformal class. Since in odd dimensions there is 
no conformal anomaly as we shall demonstrate momentarily, such a two stage 
process cannot be used. The situation can be improved considerably if we do 
not rely on the use of the two stage process just described. We would like to 
explain this fact in some detail now. 

Even though the transformational properties (with respect to conformal 
transformations) of the Laplacian, Eq.(2.4), in dimensions higher than 2 are 
rather unpleasant, fortunately, they can be considerably improved if, instead of 
the usual Laplacian, one uses the conformal (Yamabe) Laplacian D g defined by 

□ s = A g + aR(g), (6.1) 

where a — a^ 1 = jgE§- By construction, in 2 dimensions it becomes the 
usual Laplacian. In higher dimensions its transformational properties are much 
simpler than those for the usual Laplacian (e.g see Eq.(2.4)). Indeed, it can be 
shown [48] that 

n e2fg = e-^n g e^. (6.2) 

This result can be easily understood if we use the results of Section 2. Indeed, 
since e 2 ? = ip p ~ 2 and since p = -M^, we obtain at once e^ -1 )/ = ip y while the 

factor e~^2+ 1 )^ is transformed into i^ 1_p . From here, it is clear that Eq.(2.12) 
for scalar curvature can be equivalently rewritten as 

R(g) = aupi-'iAgV + a^Rig)?), (6.3) 

where g = e 2 ^ g , so that 

R{g) - an e 2f g . (6.4) 
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In practical applications it could be more useful to consider two successive 
conformal transformations made with the conformal factors e 2 t and e 2h . If 
e 2f = ip p ~ 2 and e 2h = ip p ~ 2 then, we obtain, 

A § V + or^RisW = ^ p [\ (W) + a^Rig) (W)]- (6.2b) 

This result is, of course, equivalent to Eq.(6.2a). 
Consider now the following path integral 

exp {-Tig)) = J D[<p] exp{-S Ds {<p)} (6.5) 

where 

= af d^Hv^f + a-'Rig)^} 

JM 

= a [ d*xJg[tpO g tp] = [ d d x^g~R{g) (6.6) 

JM JM 

The conformal factor ip~ v in Eq.(6.3) is eliminated by the corresponding factor 
coming from the volume factor of yfg — ^/exp(2 f)g — ip p y/g 12 . 

Corollary 6.1. Thus, Eq.(6.5) is the exact higher dimensional analog of the 
two dimensional path integral, Eq.(5.5), and, hence, problems related to higher 
dimensional CFT are those of Ricmannian 13 quantum gravity [49] . In this paper 
we are not considering the pseudo-Riemannian case relevant for true Einsteinian 
gravity. It will be discussed in a companion publication [32]. 

The question immediately arises: If the path integral, Eq.(6.5), is such an 
analog, is there a higher dimensional analog of Eq.(5.6)? The answer is "y es '\ 
if the dimension of space is even, and " no" if the dimension of space is odd [48] . 



6.2 Lack of conformal anomaly in odd dimensions 

In view of its importance, we would like to provide a sketch of the arguments 
leading to the answer "no" in dimension 3. Clearly, the same kind of arguments 
will apply in any other odd dimension. In doing so, although we follow the 
arguments of Refs. [48,50], some of our derivations are original. We begin by 
assuming that there is a one parameter family of metrics: g(x) —exp(2xf)g. 
Next, we define the operator 8j via 

S f D 9 = fa U=0 O cxp{2 xf)g ■ (6.7) 

Taking into account Eq.(6.2) we obtain explicitly 

SfD g = -2fD g (6.8) 

12 This result should be understood as follows. We have §ij = tp p ~ 2 gij -.Prom here, we have 

for dctcrminantsv's = [ i P p ^ 2 ]^ \/~§ = <P p ^/g in accord with Eq.(3.5). 
13 Sometimes called Euclidean quantum gravity. 
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and, accordingly, 
These results allow us to write for the zeta function 14 



8 f e- tu °=-t{8 f n g )e- in : (6.9) 



oo 



1 



T(s) 

-2s 



o 



dtt s Tr(-2fn g e- ta *) 



dtt s - L Tr(fe- ta °). (6.10) 



o 



The last line was obtained by performing integration by parts. Since for small 
t's it is known that, provided that dimkerD g = 15 , 

Tr(/e-* D «) ~ jr I f f(x)u k (x)dvol I t k - d '\ (6.11) 
where dvol — y/gd d x, as usual. Using this result in Eq.(6.10) produces 

1 oo 

*/Cn 9 (0) = ' dtt^Tr{fe-^)+ j ' dtt'-'Trife-^)) 

1 

OO 



Since for s — > + we have (l/r(s)) ~ s, the second (regular) term in brackets 
will become zero when multiplied by the combination j=^|y , while the first term 
will become zero even if it might acquire a pole (in case if d = 2k). Hence, for 
all dimensions d > 3 we obtain <S/Crj 9 (0) = or 

Cn B (0) = Cn a (0). (6.13) 

This result can be used further. Indeed, if we write 

*/[r(*)Cn»l =ns)[5f(nM + s5f(' n (0) + O(s 2 )} (6.14) 



14 E.g. see Eq.(5.8) and take into account the identity x~ s F(s,x) = J dyy s ~ x cxp(— xy). 



15 In the case of path integrals this is always assumed since zero modes of the corresponding 
operators are associated with some kind of translational, rotational, etc. symmetry. To elim- 
inate the undesirable dilatational symmetry, one actually should use the Yamabe functional, 
Eq.(3.6), as explained in Section 3. Although this is silently assumed thus far, arguments 
additional to those in Section 3 will be introduced further below. 
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and take into account Eq.(6.13) and the fact that sT(s) = 1 we obtain, 

oo 

*/C'n 9 (0) = S f J dtt- 1 Tr(e- tD «) | s=0 
o 

= - 2 <P 1+^1% + J dttS - lT ^- tDe )) l-o • (6-15) 

Applying the same arguments to Eq.(6.15) as those that were used for Eq.(6.12) 
we conclude that, provided dimkerD 9 = 0, in odd dimensions, <5/Ca g (0) = 0, 
that is 

Cn s (0) = CU(°)- ( 6 - 16 ) 
In view of the Eq.(5.9) this leads also to 

dctD 3 = detn § , (6.17) 

QED. 

6.3 3d CFT path integrals 
6.3.1 General remarks 

Our previously obtained results can be further refined if we recall Eq.s(5.9)- 

(5.11). In particular, let e 2 ^ in Eq.(5.11) be rewritten as some nonnegative 
constant I . Then we obtain 

Cg(s) = l s ( g (s). (6.18) 

This result is consistent with Eq.(6.13) for s = 0. Differentiation with respect 
to s produces 

^(0) = C 9 (0)lnJ + C;(0). (6.19) 

In view of Eq.(5.9), this result is equivalent to lndet D g = lndet D g — C 3 (0) Inl- 
and apparently contradicts Eq.(6.17). But the contradiction is only apparent 
in view of the earlier footnote. The situation is easily correctable if in the path 
integral, Eq.(6.5), we replace the action functional S\j(<p) by that of Yamabe 
given by Eq.s(3.2) (or (3.6)). This, by the way, allows us also to fix the value 
of C 3 (0): provided that we identify the constant I with the volume V, the value 
of C g (0) = |. After this, the situation in the present case becomes very similar 
to that encountered in the previous section. That is, instead of the functional 
T{tp) given by Eq.(5.15), we can consider now the related functional given by 

F(<p) = lndctDg -C 9 (0)lnV 

= lndetCL- -InV. (6.20) 
V 
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For the path integral calculations, a functional defined in such a way is not 
yet sufficient. To repair this deficiency we have to impose a volume constraint. 
That is we need to consider the path integral of the type 

Z Y (V) = f D[<p]6{ [ d d x^ - V)exp(-S[<p}) (6.21) 
J Jm 

with S[tp] given by Eq.(3.2) (or (3.6)). 

Remark 6.2. The path integral Zy(V) (Y in honor of Yamabe) is the 
exact higher dimensional analog of the "stringy" path integral Sl(A) given by 
Eq.(5.33). In view of Eq.(3.6), it also can be viewed as the path integral for 
pure (Euclidean) gravity in the presence of the cosmological constant 16 . 

In complete analogy with Eq.(5.32), the standard path integral for the self 
interacting scalar ip 4 (or LGW) field theory is obtainable now as follows 

oo 

J D [4>] exp{-S LGW m = J dVe- bv Z Y (V). (6.22) 

o 

But, since the variation of the Yamabe functional produces the same Eq.(3.3) 
as can be obtained with help of the L-G functional, Sl-g(4>), (which, in view 
of the Remark 3.2., can be identified with S[ip] defined by Eq.(3.7)) one can 
develop things differently but, surely, equivalently. 

To this purpose, instead of the functional S[<p] given in Eq.(3.2) we use 

Sv[<p] = A[ dVM(v^) 2 + R(g)v 2 } (6-23) 

V v JM 

and replace S[ip] in the exponent of the path integral in Eq.(6.21) by SV[y] from 
Eq.(6.23). Then, instead of Eq. (6.22), we obtain, 

OO 

J £>[0exp{-S L _ G (0)}= J dVZ Y (V), (6.24) 

o 

where the sign = means "supposedly". This is so, because, in view of the 
Remark 3.3., at the level of saddle point calculations the left hand side and 
the right hand side produce the same G-L equation. Beyond the saddle point, 
calculations are not necessarily the same. Although we plan to discuss this 
issue in subsequent publications, some special cases are further discussed below 
in this section. 



16 This fact can be used (with some caution in view of Remark 4.3.) as alternative formu- 
lation of the quantum gravity problem. 

17 To be in accord with standard texts on statistical mechanics, e.g. see [17], the constant 
scalar curvature R(g) should be identified with the squared "mass", i.e. with m 2 = a \T — T c \ 
where T and T c are respectively the temperature and critical temperature. Also the sign "-" 
in front of A term should be replaced by "+" 
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Remark 6.3. It should be clear, that at the level of the saddle point calcu- 
lation, replacement of the functional S[tp] in Eq.(6.21) by F{y>) from Eq.(6.20) 
is completely adequate, so that the sequence of steps in the analysis performed 
for the two dimensional case in Section 5 are transferable to higher dimensions 
without change. 

Remark 6.4. Although in 3 dimensions we have the result given by Eq.(6.17), 
which forbids the use of identities like that in Eq.(5.6) still, based on arguments 
just presented, the functional !F(ip) defined by Eq.(6.20) should be used in the 
exponent of the corresponding path integral replacing that given in Eq.(5.15) in 
2 dimensions. Since by doing so one will be confronted with the same type of 
minimization problems as discussed earlier in Section 5 18 , the functional integral 
thus defined is an exact 3 dimensional analog of the path integral, Eq.(5.30). 



6.3.2 A sketch of 3d CFT calculations 

Our obtained results allow us to proceed in complete accord with Section 5. 
Thus, let g(u) = e^g be a 1-paramctcr family of metrics of fixed volume 
and such that g(0) = g. This implies that 0(0) = and / e^dVo = V. 

M 

In two dimensions, using results of OPS [41], especially their Eq.(1.12), it is 
straightforward to obtain the following result 

|-(- lndet 'A fl(u) ) U =0 = C>) = / j>K(g)dV g (6.25) 

M 

where K(g) is the Gaussian curvature for the metric g. C 9 (0) and represent 
SrCg(tt)(0) U=o an d -§^4>{ u ) |«=o r espectively. Volume conservation implies 

^ je+WdV a |„=o= j 4>dV g = (6.26) 

M M 

in accord with the earlier result, Eq.(5.17). If the Gaussian curvature K(g) is 
constant, then Eq.(6.25) and Eq.(6.26) produce the same result. This implies 
that J^C s (tt)(0) | M =o= 0, that is g is the "critical" (extremal) metric. In view 
of Eq.(6.25) this also means that, for such a metric, the free energy (e.g. see 
definition given by Eq.(5.5)) attains its extremum. By combining Eq.s (6.10), 
(6.14) and (6.15) we obtain at once 

oo 

*/Cn,(0) = S f J dtt^Trie-^) | a=0 



oc 

= J dtt s Tr(-2fn g e- ta *)\ s=0 . (6.27) 
o 

18 This conclusion was reached without any reliance on path integrals and on physical ap- 
plications in Ref.[51]. In this and related Rcf. [20] the extremal properties of determinants of 
□ g and A 9 with respect to variations of the background metric were studied. 



24 



From here we obtain essentially the same result as the main theorem by Richard- 
son [20], i.e. 

C' n . (0) \ u=0 = = J dV g j>(x)n g ((l, x, x). (6.28) 

M 

That is, provided that we replace D g by A g and require that locally ((1, x, x) = 
const, with the heat kernel Q(s,x,x) given (as usual) by 

C{S , X , X ) = J2^, (6.29) 

k=i A fe 

we obtain the main result by Richardson, Ref.[20], e.g. see his Theorem 1 
and Corollary 1.1. Here ^Pk( x ) are eigenfunctions of the Laplacian (or Yamabe 
Laplacian respectively) corresponding to eigenvalues \k with = Ao < Ai < 
A2 < • • • < Afc < • • •. In the two dimensional case, the condition for critical- 
ity: C g (0) = is local meaning that, provided the volume is constrained, the 
constancy of the Gaussian curvature K (g) at a given point of M is caused 
by the metric g for which C s (0) is extremal. In 3 and higher dimensions, con- 
stancy of the curvature at the point of M is being replaced by constancy of 
C(l,x,x) under the same conditions of volume conservation. Such a condition 
is only necessary but is now not sufficient, since the analog of Moser-Trudingcr 
inequality, Eq.(5.21), used to prove sufficiency in 2 dimensions does not exist. 
Instead, one should study locally the second variation of C§(0) with respect to 
the underlying background metric in order to decide if such a (local) extremum 
is a maximum or minimum. Fortunately, this task was accomplished in Ref.s[20, 
51]. In particular, Richardson [20] obtained the following theorem 

Theorem 6.5. (Richardson) The Euclidean metric on a cubic S-torus is a 
local maximum of determinant of the Laplacian with respect to fixed-volume 
conformal variations of the metric. 

Remark 6.6. This Theorem is proven only for the cubic 3-torus. The 
word "local" means that there could be (or, there are, as we shall demonstrate) 
other 3-tori also providing local maxima for determinants. In fact, according 
to the result by Chiu [52], all determinants of flat 3-tori possess local maxima 
so that the determinant for the face centered cubic (fee) lattice has the largest 
determinant. Physical implications of this fact can be found in our recent work, 
Ref.[47], and also will be briefly considered in Section 8 below. 

Remark 6.7. By using the same arguments as in Hawking's paper, Ref. 
[53], the above theorem by Richardson can be given a physical interpretation. 
It is based on use of the saddle point methods applied to the Yamabe path 
integral, Eq.s (6.21), (6.24). 
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The second variation of the Yamabe functional was calculated by Muto [54] 
(see also Ref.[49]) with the result: 

=^[Jdv g (a(v^) 2 -R(g)^)l (6-30) 

where the constant a = d — 1. As in the case of quadratic actions in the flat 
space [17, 36] the second variation (with volume constrained to be equal to 
one) looks very much the same as the original quadratic Yamabe functional, 
except for the "wrong" sign in front of the scalar curvature. This "wrong" 
sign has important physical significance in both statistical mechanics and high 
energy physics. While the significance to high energy physics will be discussed 
in the companion publication, Ref.[32], here we would like to comment on the 
significance to statistical mechanics. For this we need to recall that for the G-L 
free energy functional F{<p} given by [17] 

F{<p} =K + \j d*{c (V^(x)) 2 + V(x) + ^ 4 (x)} (6.31) 
the second variation is given by 



6 



2 



(6.32) 



_5<p(x)6<p(-x') h=o 
where 

Mx)Mx') ^ M lv=v " = ( ~ c y2 +a + 3blfl) S{x ~ x °- (6 ' 33) 

with ip being obtained from the extremum condition 

5 -Pr^ = -c\/ 2 ^ + a^ + bp 3 -h = Q (6.34) 

In a special case of a coordinate independent field ip Eq.(6.34) is reduced to 

cup + bipl = 0. (6.35) 

For a < we obtain : cp = ±(|a| /6)5 and ip = 0. The solution <p = is not a 
minimum for the free energy and, hence, should be discarded. Of the two other 
solutions, the system chooses one of them (which is interpreted in literature as 
spontaneous symmetry breaking). They both have the same free energy 19 . If we 
let now y>(x) — ip Q + 5f = ip + jj(x) and Taylor series expand ^F{f} the result 
will coincide with Eq.s(6.32) and (6.33). Notice that for a < a combination 
-|o| + 36(^g = 2 \a\ = m 2 (as compared to a in the high temperature phase). 
So, if in the high temperature phase to 2 = R(g) with R(g) > 0, in the low 



19 Once the choice is made, the order parameter can be considered in all subsequent calcu- 
lations as positive. 
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temperature phase (i.e. below the criticality) R(g) < 0! Moreover, by choosing 
d = 4 in Eq.(6.30) we obtain a combination: m 2 = -^-R(g) = 2R(g) in accord 
with the standard G-L results. The constant c in Eq.(6.33) is left unspecified in 
G-L theory while now it is equal to 6 (for d = 4). The spaces of constant negative 
curvature are hyperbolic. 

Remark 6.8. In our earlier work, Ref[55] we emphasized the major im- 
portance of hypcrbolicity in statistical mechanical calculations related to the 
AdS/CFT correspondence. In this work we arrived at the same conclusions 
quite independently. 

From the work by Muto [54] we know that: a) if R{g) is positive, the sec- 
ond variation can be made positive for appropriately chosen ip 20 , b) if R(g) is 
negative, the second variation is positive for the same reasons. The positivity 
of the second variation implies that the extremal constant curvature metric g 
provides a locally stable minimum for 1Z(g(t)) defined by Eq.(4.9). That is, 
using results of Section 4, the Einstein metric obtained as solution to Eq.(4.2) 
is stable among nearby metrics. This conclusion will be reinforced in the next 
(dynamical) section. 

Remark 6.9. It is interesting to notice that the calculation of higher or- 
der fluctuation corrections to the Yamabe path integral, Eq.(6.21), involves 
calculations on the moduli space of Einsteinian metrics, Ref.[56]. This obser- 
vation provides a strong link between higher dimensional LGW theory and two 
dimensional string inspired CFT's discussed in the previous section. Natu- 
rally, Eq.(6.24) can be used to investigate to what extent the final results of 
conventional field-theoretic calculations, Ref. [17,36], might differ from more so- 
phisticated string-theoretic calculations in the style of Ref.s[14, 15, 46]. Additional 
details can also be found in our work, Ref. [47]. 

7 Critical dynamics and the Yamabe and Ricci 
flows 

7.1 Physical motivation 

The results obtained in previous sections are somewhat unrealistic from a purely 
practical (physical) point of view. The situation in our case is analogous to that 
known in thermodynamics. Recall, that this discipline emerged from practical 
needs to improve the efficiency of heat engines. Clearly, rigorously speaking, it 
is not applicable to such devices since, by definition, it is valid only for time- 
independent phenomena. Hence, it cannot be used at all because it usually 
takes a very long (if not infinite) time for the system to equilibrate, especially 
near the critical temperature T c . Thus, the description of phase transitions by 
the G-L theory is valid only in a rather narrow range of temperatures around 

20 This can be easily understood if we expand ip into Fourier series made of eigenfunctions of 
the Laplacian and take into account that for any closed Ricmannian manifold the spectrum 
of the Laplacian is nonncgative and nondecreasing [29]. 
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T c . To bring a physical system into this range of temperatures (under constant 
pressure) requires varying the temperature of the surrounding environment in 
time. In addition, by definition, the "true" phase transition should take place 
only in the thermodynamic limit (of infinite volume with particle density kept 
constant). Since in the real world the systems under study are always of fi- 
nite size, this requirement is implemented by imposing physically appropriate 
boundary conditions, e.g. periodic. Use of some appropriate boundary condi- 
tions causes systems undergoing phase transitions to actually "live" on some 
manifolds/orbifolds. Under such circumstances the topology and physics be- 
come intertwined. The signature of such boundary effects can be seen already 
in the calculation of determinants discussed in Sections 5 and 6. More on this 
subject is discussed in our earlier work, Ref.[55], where, in accord with Re- 
mark 6.8., we emphasized the role of hyperbolic spaces in the theory of phase 
transitions 21 . 

The equilibration process known in the physics literature as critical dynamics 
can take a very long time (infinite in the thermodynamic limit for T => T c ). 
Since idealization of reality is typical in physics, in this work we adopt the 
pragmatic (physical) point of view by considering systems of finite size with 
appropriately chosen boundary conditions. Then, the results of the previous 
sections are limiting cases of more general time-dependent G-L theory consid- 
ered for such systems. 

Development of time-dependent G-L theory was initiated by Landau and 
Khalatnikov in 1954 [57] and is also phenomenological. As such, it is based on 
the assumption that an order parameter ip satisfies the relaxation equation of 
the type given by 

with functional F{>p} defined by Eq.(6.31) and the "friction" coefficient 7 is 
some (assumed to be) known non-negative constant. By rescaling time it can 
be eliminated. Such rescaling is assumed in all calculations below. Since such 
an equation was postulated, its validity was checked by real experiments with 
excellent outcome [58]. Being armed with such results, we would like to develop 
the mathematical formalism of the previous sections to account for the effects 
of critical dynamics. 

7.2 Mathematical motivation 

Eq.(7.1) is an example of the gradient flow. In this subsection we would like to 
place the Landau-Khalatnikov theory of critical dynamics on a more rigorous 
mathematical basis also involving the notion of the gradient flow. 

We begin with reviewing the results of OPS, Rcf. [41], (obtained for two 
dimensional manifolds without boundaries) in the style of Hamilton's work, 

21 Remark 6.8. is not in contradiction with the remark 6.6. This is explained in great detail 
in our recent work, Ref.[47] and, is also discussed in Section 8. 



28 



Ref.[16]. To this purpose we notice that Eq.(5.19) is an extremum of the func- 
tional F(ip), Eq.(5.14). Clearly, if we write ip = ip + eh in this functional so 
that 

/ dAgh exp(2ip) 

^+eh) \^=-fdA gM)h+ \ JdA g R( 9 )h-2, X (M r jdAgexpm 

MM M 

(7.2) 

then, following OSP, the scalar product (, ) in the tangent space to each point 
at which ip is defined can be introduced according to convention: 

{gradT, h) = {5T) h = ^-T{^ + eh) | e=0 . (7.3) 

In accord with Eq.(3.6) by OPS, Ref.[41], this convention implies that 

1 2^ X (M)exp(2^) 

2 {9) /exp(2^)dA fl ' 



M 

This result is to be compared with Eq.(5.19). Naturally, Eq.(5.19) corresponds 
to the condition of equilibrium, i.e. gradT {ip} = 0. Away from equilibrium 
the dynamical equation postulated by OPS (see below) coincides with Eq.(7.1) 
postulated much earlier by Landau and Khalatnikov. To see this, we need to 
take into account Eq.s (5.2), (5.19) and (5. 20). With their help we would like to 
rewrite Eq.(7.4) in the following equivalent form, 

-gradf^} =- ( < R(g) > -R(g))e 2 ^ , (7.5) 

where < R(g) > 



( j dAR(g) ) / J dA with dk given by Eq.(5.3). Accord- 

\M / M 



ingly, in view of Eq.(7.1), up to time rescaling, the equation describing critical 
dynamics can be written as follows 

^ =< R(g) > -R(g) (7.6) 
coinciding with Eq.(3.9) of OPS 22 

Suppose that the equilibrium (the fixed point) condition, 

< R(g) > -R(g) = 0, (7.7) 

is valid in dimensions higher than two as well. Then it coincides formally with 
Eq.(2.12) defining scalar curvature or with Eq.(4.12) where it was obtained 
from the Hilbert-Einstein functional for pure gravity in the presence of the 
cosmological term. In view of this observation, we would like to demonstrate 
that equations like Eq.(7.6) whose fixed points are determined by the equations 
of the type given by Eq.(7.7) indeed can be obtained in higher dimensions. In 
this work only 3 dimensional results will be discussed in some detail. 

22 In making such a comparison (as in Section 5) we took into account that R = 2K. 
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7.3 Critical dynamics in dimensions higher than two: 
Yamabe vs Ricci flows 

Following Hamilton [16], we begin with 

Definition 7.1. The normalized Ricci flow is described by the dynamical 
equation given by 

d 2 

■EZ9H = ~ A m < R(9) > -2Rij{g). (7.8) 



In writing this equation we use the notations of Section 4 where we defined 

< R(g) >= (Im R Va ddx ) / Im V9 ddx = (S M Mn) I f M d »- The above now 

equation formally exists for d > 2 and should be considered along with some 
prescribed initial condition: gij(t = 0) = §ij. In two dimensions it is always pos- 
sible to write Rij(g) — \R-9ij so that all two dimensional spaces are Einsteinian 
[16]. Because of this fact, Eq.(7.8) is converted into 

^ gij - (< R(g) > -R(g)) gij , d = 2. (7.9) 

Definition 7.2. Eq.(7.9) defines the Yamabe flow. As such, it can be 
considered for any d > 2. 

Remark 7.3. If, like in Section 2, we choose g^ = e^g^ (0) and substitute 
this result into Eq.(7.9), we reobtain Eq.(7.6). Hence, Eq.(7.6) does describe 
the Yamabe flow. 

Corollary 7.4. Using results of previous subsection (at least in two dimen- 
sions) the Yamabe flow is the gradient flow. 

Corollary 7. 5. Since the Yamabe flow in Eq.(7.9) is defined for any d > 2, 
the fixed points for such a flow should coincide with those given by Eq.(7.7). 

Remark 7.6. Only in 2 dimensions are the Yamabe and the normalized 
Ricci flows (essentially) equivalent. This is explained in great detail by Hamilton 
[16]. Since results of this paper crucially depend on the observation that only in 
2 dimensions it is always possible to write Rij{g) = \Rgij, it should be clear that 
attempts to extend this equivalence to higher dimensions are destined for failure 
in those cases where the flow takes place on spaces which are not of Einstein 
type. Hence, in this case the task lies in finding conditions under which the Ricci 
flow, for which the initial metric is chosen to be of not necessarily Einstein-type, 
leads to an Einstein-type metric(s) as the fixed point(s) solution(s) for such a 
flow. This is discussed in the next subsection. 

Even though Yamabe and Ricci flows are different in higher dimensions, one 
may still ask a related question : can one use some Ricci flows in order to 
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obtain results for Yamabe -type flows? We would like to demonstrate that this 
is indeed possible. To this purpose, using Eq.(7.9) we obtain, 

\g ii ^9i j =<R{g)>-R{9). (7.10) 

In addition, since in any dimension -§^dfi = ^/det gijd d x — \g l i -^g^dfj, , we 
obtain as well 

§1 Im d ^ = Im ^ = J M (< R > -R{g)W = 0. (7.11) 

This equation is compatible with Eq.(7.8) and valid for d > 2. Indeed, if we 
take into account that -§^d^i — \g l i -^g^dfj,, use Eq.(7.8) for ^gij, multiply both 
sides of this equation by g^dfi and integrate the obtained result over M, we 
reobtain Eq.(7.11). Hence, after some manipulations Eq.(7.8) describing Ricci 
flow (and valid for d > 2 ) can be brought into a form identical with Eq.(7.10) 
describing Yamabe flow (obtained originally for d = 2). This observation does 
not imply that both flows are equivalent as Remark 7.6. indicates. Nevertheless, 
the Yamabe flow can be looked upon as some special case of the Ricci flow. 

Corolary 7.7. The arguments just presented demonstrate that solutions for 
the Yamabe flow can be obtained from solutions for the normalized Ricci flow. 
Moreover, in both cases the fixed points for such flows are given by Eq.(7.7) 
thus providing justification for the compatibility of 2 and 3 dimensional CFT 
discussed in previous sections and, hence, for the unique pathway for extending 
the results of 2d CFT to 3 and higher dimensions. 

Corollary 7.8. Since the fixed point, Eq.(7.7), is the traced form of the 
fixed point equation for the Ricci flow, Eq.(7.8), that is of the equation 

^ 9ij <R(g)>-R ij (g)=0, (7.12) 

which is equivalent to Eq.(4.2), this means: a) that the fixed point solutions 
of the normalized Ricci flow always produce the Einstein spaces, provided that 
these fixed points are stable and, b) that the arguments presented in Remark 
4.3. are indeed valid. 

7.4 Stability of the fixed points solutions 

The results obtained previously in this section cannot be used until the stability 
analysis of these solutions is performed. Since according to the Corollary 7.5. 
solutions for the Yamabe flow can be obtained from those for the Ricci flow, it is 
sufficient, in principle, to consider only the stability of the fixed point solutions 
for the Ricci flows. This path is not the most physically illuminating however 
as we would like to explain now. 

7.4.1 Dynamics of the Yamabe flows 

The Yamabe flow described by Eq.(7.9) should be supplemented by the ini- 
tial condition, e.g. gij(t = 0) = gij(0) = g. In order to study the evolution 
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of this metric it is convenient to write it in the following form [59]: gij{t) = 
[ip(t)}^ g {j (0). Using such a substitution in Eq.(7.9) (along with Eq.(2.12)) al- 
lows us to rewrite it in the following equivalent form 

^ - -^ 2 ">A^ + R{g) V )+ < R(g) > V (7.13) 

to be compared with the phenomenological result, Eq(7.1), by Landau and Kha- 
latnikov. Although these equations look somewhat different, their fixed points 
(if these are stable) are the same 23 . In view of the results of previous sections, 
we believe that the equations of critical dynamics known in the physical lit- 
erature, e.g. see Ref. [58], all having their origin in the work by Landau and 
Khalatnikov [57], should be replaced by Eq.(7.13) also describing critical dy- 
namics of physical systems of the Ginzburg-Landau type. Such a replacement is 
in accord with the requirement of conformal invariance at criticality as described 
in detail in our earlier work on AdS/CFT correspondence [55]. 

Remark 7.9. Evidently, physical realization of the Yamabe and Ricci 
flows 24 depends crucially on the possibility of such a replacement and, in view 
of results of this work, it can be considered as proven. 

Using Eq.(7.10) we would like to introduce the following quantity 

V = ^(< R(g) > -R(g))- (7-14) 

In view of Eq.(7.11), integrating Eq.(7.14) we obtain, 

f M Vd t i = 0. (7.15) 

By combining Eq.s(7.11) and (7.13), after some calculations which can be found 
in Ref. [59], we obtain, 

| < R(g(t)) >= ^ J M R(g) V d». (7.16) 
Next, by combining Eq.s(7.14) and (7.16) we further obtain, 

^ < R(g(t))) >=2f M R(g)(< R(g) > -R(g)) = -2 / M (< R(g) > -Rig)) 2 , 

(7.17) 

where the second equality comes from the constraint, Eq.(7.15). The last ex- 
pression possesses an entropic meaning. Indeed, following Boltzmann [60], we 
introduce the if-function (the entropy) for some gas made of hard spheres elas- 
tically scattering from each other . Explicitly, it is given by 

H = - J flnfdv, (7.18) 

23 In two and three dimensions these fixed points were discussed in previous sections. 
24 And, ultimately, the physically assisted proof of the Poincare' and geometrization conjec- 
tures depend crucially on such an assumption. 



32 



where the distribution function /(v,i) obeys the Bolzmann's equation 

§- t f = JW^/K,*) - /(v,i)/( Vl ,i)) |vi - v| J((9, |vi - v|)dvidn (7.19) 

with |vi — v|) being some known function, dQ — sa\9d9d<p, and v' = 
v + n(n • g),Vj = vi — n(n • g) being respective velocities of the colliding 
particles after (primed) and before (nonprimcd) scattering on each other, here 
g = Vi — v, and the vector n is a unit normal 25 . Based on these results, it 
can be shown that dvdvi — dv'dv^. Using these results, let us consider cal- 
culation of JjiJ.By combining Eq.s(7.18) and (7.19) one obtains after some 
calculation -j^H > with an equality only for the equilibrium case for which 
f*(-v',t)f*(v[,t)-f*(v,t)f*(-v 1 ,t) = 0. This condition is leading to the Maxwell's 
distribution of particle velocities, which produce a well tested equation of state 
for the ideal gas. In accord with thermodynamics, such defined entropy riches 
its maximum at equilibrium. 

Remark 7.10. The obtained equilibrium result for Maxwell's distribution 
depends crucially on the assumption of occurrence of only binary collisions in 
the gas. Boltzmann's equation, Eq.(7.19), reflects just this fact. The obtained 
results become much less tractable if one would like to account for ternary 
and higher order collisions. The same is true if one uses the full renormaliza- 
tion group flow equations [4] instead of just the leading order terms in these 
equations. The proof of the Bolzmann's H-theorem as well as the proof of the 
Poincare' and geometrization conjectures are based on the validity of such ap- 
proximations. The validity is assured by excellent agreement with experiment 
in Boltzmann's case and with equally good agreement with experiment in the 
case of Landau-Ginzburg critical dynamics described by the Yamabe-type flow. 

Based on these physical arguments, we may choose — < R(g(t))) > as our 
entropy (provided that the volume of the system is normalized, say, to unity). 
Then, we can identify Eq.(7.9) (or (7.13)) with a Boltzmann-type equation. Our 
previously obtained result, Eq.(4.12), then provides the physically meaningful 
equilibrium solution. Moreover, using these results we arrive at the conclusion 
that the fixed point solutions for the Yamabe flow are described by the Einstein 
equation, Eq.(4.2), leading to the Einstein manifolds of constant scalar curvature 
in accord with the Remark 4.3. 

7.4.2 Dynamics of the Ricci flows 

In this subsection we would like to demonstrate in some detail that, indeed, 
some Ricci flows produce results that are in agreement with those obtained 
for the Yamabe flows. For this purpose, following the logic of the previous 
subsection we have to find an analog of the entropy for the Ricci flows. This 
task was accomplished by Perelman in the first of his 3 papers, Ref.[l]. He 
noticed, albeit indirectly, that, unlike the Yamabe flow (which is the gradient- 
type flow), the Ricci flow is not gradient. The proof can be found in Ref.[61]. 

25 See [60] for details. 
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Accordingly, one cannot find an entropy for such a flow and, hence, to study 
its stability. Pcrclman found ingenious way around these difficulties. In view of 
the existing literature [2,3, 10, 61-63], our exposition will be brief emphasizing 
the physical aspects of the results he had obtained 26 . 

The Boltzmann-Nash entropy and the dilaton gravity. In the case of 
the Yamabe flow our choice of — < R(g(t))) > as an entropy seemed somewhat 
artificial in the sense that it was not motivated by any systematic procedure 
involving entropy as such. We would like to correct this deficiency in this sub- 
section. We begin with the following 

Definition 7. 11. The heat operator □ := d t — A g has an operator □* = 
— dt — A g + R(g) as its adjoint with respect to the Ricci flow. E.g. see Ref.[61], 
page 22. 

Remark 7.12. The proof is based on the observation that upon some 
rescaling [5] it is possible to supress the normalization in the Ricci flow, Eq.(7.8), 
that is to eliminate the factor ^9ij < R{g) >■ As result of this elimination, 
Eq.(7.10) acquires the form 

\j i ^9a = -Ris) (7.20) 

used in the proof of the conjugacy of operators □ and □*. 

Let low gij be a solution of the non normalized Ricci flow equation and 
u = exp(-f) be the solution of the equation D*u = for some function / 
such that . J M e~f dV = 1. By nanalogy with Boltzmann's entropy, Eq.(7.18), 
following Rcf.[61] we define the Nash entropy N(u) as follows: 

N(u) =J M ulnudV. (7.21) 

The time derivative d t N(u) can be calculated now with help of the adjoint 
equation D*u = 0. Indeed, we obtain: 

d t N(u) = j M {d t u\nu dV + d t u dV +ulnu d t dV) 
= J M ((d t - R(g))ulnu + d t u )dV 
= f M (-A g ulnu+{R(g)-A g )u)dV. (7.22a) 

On a closed manifold, the integral of A g u vanishes. Because of this, the last 
result acquires the form 

d t N(u) = J M (^L + R{g)u)dV = J M (\v g f\ + R(g))e- f dV = F( 9ij ,f) 

(7.22b) 

The functional T{gij, /) is Perelman's entropy [1] and is also an action for the 
dilaton gravity [ 8,64]. When / is a constant, taking into account the sign 
differences between the Bolzmann's entropy, Eq.(7.18), and the Nash entropy, 

26 Naturally, only small part of Perelman's results is discussed in this paper. 
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Eq.(7.21), we obtain back our earlier result for the entropy — < R(g(t)j) > as 
required. 

Remark 7.13. The entropic nature of the Einstein- Hilbcrt and dilaton 
gravity actions just described might be especially useful for applications to black 
hole thermodynamics. Indeed, the attempts to use the Ricci flow for the de- 
scription of black hole dynamics were recently made in Ref.[65]. We plan to 
return to this issue in subsequent publications. 

Remark 7.14. If one starts with the functional ^{gij.j) without prior 
remarks on its entropic origin, one might employ the dilaton field / which is 
time- independent. This leads to some technical simplifications to be discussed 
below. 

Having defined the entropy for the Ricci flow, the task now lies in obtaining 
the analog of Eq.(7.17). We would like to do so in a way consistent with our 
results for the static case. 



Steady Solitons. In particular, following the logic of Section 7.2., we should 
begin with the replacement of the two dimensional functional T{ip) by its three 
dimensional analog, e.g. by S[ip] defined by Eq.(3.7), or, equivalently, by S[ip] 
defined by Eq.(3.2). To make a connection with Perclman's work, we need 
temporarily to make a short but important detour. To this purpose we set 
(p = e~$l 2 in Eq.(7.22). This causes the functional ^F(gij,f) to acquire the 
following new look 

<p) = Jm( 4 + R(9)V 2 W (7.23) 

to be compared with the functional E[tp] defined in Eq.(3.2). By analogy with 
Eq.(3.4), now we can define the constant X g in terms of the Raleigh quotient, 
i.e. 

T J M ^dV 

In accord with Eq.(3.3), the constant X g serves as an eigenvalue in the equation 27 



A ff = inf^^. (7.24) 



AAglp + R{g)y - \ g (p = 0, (7.25) 

where ip is the minimizer for the Raleigh quotient, Eq.(7.24). Equivalently, 

A, = inf{/ M (4 |V^| 2 + R(g)<p 2 )dV, J M <p 2 dV = 1}. (7.26) 

By analogy with Eq.s (4.7)-(4.10), we introduce the family of metrics gij(s) = 
g-ij + $ hij in order to consider the Raleigh quotient under such a variation. 
Then, instead of Eq.(4.10), we obtain 28 

^Gfe(s)) - Jm(-M(^ + V&jf)e-ldV. (7.27) 

27 See the second footnote for the sign convention for the Laplacian. 
28 For details of such calculations, please, consult [10] and [66]. 
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Since this result is analogous to that given in Eq.(7.16), the course of action is 
going to be the same. In particular, following Pcrclman, and also [10], we notice 
that the variation of the functional F(gij,f) leads to the following coupled 
equations for the modified Ricci flow 29 

^gij = -2(R t3 + ViVj-/), (7.28a) 

^-f = -R + Af. (4.28b) 

Definition 7.15. The flow defined by these coupled equations is called the 
generalized Ricci flow. 

Finally, after making the identification -§^gij — hij and using Eq. (7.28a) in 
Eq.(7.27) we obtain the monotonicity result of Perelman 

j- s H9ij(s)) =2J M \R l0 + V.V./I 2 e-fdV (7.29) 

to be compared with Eq.(7.17). In view of this comparison, we conclude that 
j^X(gij(s)) = only in the case when 

R ij + V i V j f = 0. (7.30) 

Definition7.16. In the existing terminology Eq.(7.30) is known as the 
equation for the gradient steady soliton [ 1-3, 66]. 

Following Perelman [part I] we can extend this result to the case of the so 
called breathers. 

Definition 7.17. A metric gij(t) evolving by the Ricci flow is called a 
breather, if for some t\ < t2 and a > the metric agij(t\) and gijfa) differ 
only by a diffeomorphism. The cases when a = 1, a < 1 and a > 1 correspond 
to steady, shrinking and expanding breathers respectively. Solitons are trivial 
breathers for which the above relationship between metrics holds for any pair 
of ti and i 2 • 

Corollary 7.18. If one considers the Ricci flow as a dynamical system on 
the space of Riemannian metrics modulo diffeomorphisms, then the breathers 
and the solitons are respectively the periodic orbits and fixed points for such a 
system. 

Corollary 7.19. In view of the result, Eq.(7.29), one can write ^(gij(si)) < 
A(5y (s2)). Because of Eq.(7.30) we obtain as well: X(gij (si)) = X(gij(s2)) Vs e 
[si, 82] - Hence, Eq.(7.30) providing a minimum for X(gij(s)) is a steady breather 
and a steady soliton at the same time. That is, there are no breathers for the 
Ricci flows in accord with Perelman [1] . 

Going back to Eq.(7.30), multiplying it by g % i and making a summation 
over repeated indices we obtain, 

R = A g f, (7.31) 

29 Please, see again the sign convention for the Laplacian in the second footnote. 
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where, as before, we used the fact that — A g f = V J Vi/. Taking into account 
arguments which lead from Eq.(7.22a) to (7.22b) we conclude that for compact 
manifolds / = const and, hence, R = 0. 

Expanding and shrinking solitons. The obtained result is a special case 
of more general result by Ivey [67] who extended the work by Hamilton [16] for 
surface to Ricci flows on compact 3- manifolds 30 . It is instructive to conclude 
this subsection by discussing the physical significance of their results, especially 
in connection with results obtained in the previous section. 

To this purpose, let us recall that the fundamental solution of the heat 
equation in d— dimensional Euclidean space can be written as [63] 

u(x, y, t) = (47tt)~^ exp(- \x - yf /4r), (7.32) 

where t = t — TorT — t = —t = (j depending upon wether we are deal- 
ing with the forward or the backward heat equation 31 . If we let y — 0, then 
such a solution can be interpreted as probability 32 since by design we have 
J M udV = 1. Using this fact, one can look for a solution of the heat equation 

_d 

on some Ricmannian manifold M by employing the ansatz u = (47rr) 2 e~' . 
This amounts of redefining the earlier introduced function /:/ = /— 5 ln(47rr) 

so that e~* = (47tt) _7 e~ s . Evidently, V/ = V/ and, therefore, A g f = A g f 
as well. To simplify matters, we would like to consider how things are done in 
the flat case. Using Eq.(7.32) we begin with the calculation of the Nash entropy. 
A simple calculation produces: Nfi at = — | — |ln(47rr). Using this result, we 
now define the properly normalized Nash entropy, i.e. N(u) — Nfi at . Explicitly, 
we obtain, 

N(u) = N(u) - N flat = J M (-f + ^)udV. (7.33) 

Perelman does not explain how he had obtained his result in his work. He 
calls such a normalized Nash entropy a "partition function", e.g. see Section 
5 of Part I of Ref.[l]. This is a bit misleading however because in the same 
section he earlier defines the partition function correctly, i.e. in accord with the 
accepted rules of statistical mechanics. According to these rules one defines the 
free energy T via \nN(u) = —j3T with [3 being subsequently identified with 
the inverse temperature (provided that the system of units is used in which the 
Boltzmann's constant, fes = 1). In the present case the role of temperature is 
played by t. Using these definitions, one can define an "energy" U =< E > in 
a familiar way as U = — ^lnTV(u) = r ^ In N(u). In our case, explicitly, we 
obtain: 

U = r 2 ^[J M ulnudV-N flat ] = J M T 2 -^(ulnudV) + ^rf M udV. (7.34a) 

30 Later, in Rcf. [68], Hamilton came up with yet another proof of Ivey's results. 
31 Here T is some pre assigned time. 

32 The fundamental solution, E.q.(7.32), also has the meaning of probability in polymer 
physics [69]. This fact is briefly discussed in Section 8. 
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33 To proceed, we need to employ Eq.s(7.22a) and (7.22b). This leads us to 
U = T 2 ^-[J M ulnudV-N flat }=J M [r 2 (R+(Vf) 2 ) + UudV. (7.34b) 



This result differs slightly from that obtained by Perclman. If we are inter- 
ested in calculation of the entropy, such a difference is useful. Indeed, since 
thermodynamically the relationship 

(3U -(3T = S (7.35) 

determines the entropy S, using already obtained results for U and T we obtain, 

S+ = J m [t(R + (V/) 2 ) - / + d]udV. (7.36) 

This result coincides with the result for the entropy of the Ricci expanders ob- 
tained in Ref.[61] where a considerably more cumbersome and lengthy pathway 
was chosen in order to obtain it. In order to obtain the entropy for the Ricci 
shrinkcrs, it is sufficient to change signs when taking time derivatives. Hence, 
we obtain at once: 

S- = f M [a(R + (V/) 2 ) + / - d]udV, (7.37) 

again in accord with Ref.[61]. At this point, one can proceed either by computing 
the heat capacity C v — {-§fU )y under the constant volume or repeat the 
arguments for the steady solitons adapted to the present case. The last pathway 
is discussed in the paper by Cao, Ref.[66]. More physically attractive, however, 
is to follow the logic of Perelman and to take into account that 34 

or 

T*C v = ^\nN{u). (7.38) 

A straightforward but lengthy calculation analogous to that used in Eq.(7.34) 
finally leads to 



r 2 a = r 4 / 



hi 



1 2 



Rij +ViVjf - — gtj 



dV (7.39) 



in accord with earlier obtained result, Eq.(7.29), for steady soliton (obtained 
from this expression in the limit |t| — > oo). The gradient shrinking (or expand- 
ing) solitons are respectively solutions to 

Rij + ViVjf - ^ gij = (7.40a) 



33 Here we took into account that f M udV = 1 as required. 
34 As is well known from statistical mechanics, 

C v =f3 2 (< E 2 > - <E > 2 ), 

where, as before, we put fcg = 1. Hence, what Perclman is simply calling a fluctuation in 
Section 5 of Part I of his work, Ref.[l], is actually a heat capacity. 
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or 

% + ViVjf + l~ gij = (7.40b) 
At 

equations. By analogy with Eq.(7.30) we multiply both of these equations by 
and sum over repeated indices in order to obtain 

R-A g f-^-=0 (7.41a) 

and 

R-b g f+± = (7.41b) 

For compact manifolds A g f = as before. As a result, in both cases we 
obtain back Eq.(4.3a). By combining these results with those which follow 
from Eq.(7.31) we thus have rederived the following result by Ivey. 

Theorem 7. 20. There are no three-dimensional solitons or breathers on a 
compact connected 3-manifold other than those of constant curvature metrics. 

Remark 7.21. This theorem was originally proved by Ivey, Ref[67], who 
was inspired by the earlier obtained result by Hamilton [16] for surfaces. Our 
derivation, however, is inspired by Perelman, Ref.[l], who uses physical argu- 
ments. Because these arguments are not presented in sufficient detail in his 
papers, they were left unappreciated in the works by other mathematicians 
[2,3,10,61-63]. 

Corollary 7.22. In view of the obtained results it appears that, at least for 
compact connected 3-manifolds, one can extract the needed physical information 
for the Yamabe How from that for the Ricci How. We refer our readers to works 
by Hamilton [5,16,70] where many illuminating additional details can be found. 
The obtained result provide needed support to claims made in Remark 7.6. 

Corollary 7.23. The generality of arguments used for reproving the Theo- 
rem 7.20. are such that they can be applied, in principle, to manifolds of any 
dimension d>2. Because of this, for compact manifolds, the Euclidean dilaton 
gravity described by the Eq. (7.22b) for any d > 2 is reduced to the more familiar 
Euclidean gravity. 



8 Discussion 

8.1 Other physical processes whose dynamics is described 
by Ricci flows 

The goal of this paper is to find some processes taking place in the real world 
providing a physical justification for Ricci flows. In the previous sections such 
processes were found, e.g. in the form of critical dynamics. The question arises 
wether there are other physical processes which also can be used for justifying 
othe existence of Ricci flows. Earlier, in the text, in Ref. [65], we mentioned the 
dynamics of black holes. In addition, we would like to notice that the fundamen- 
tal solution, Eq.(7.32), of the heat equation can be interpreted as an end-to-end 
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distanccdistribution function used in polymer physics for the computation of 
various averages, e.g. the mean square end-to- end distance for the flexible 
polymer being modelled by a random walk [69] . Under such circumstances, the 
time r is interpreted as polymer's contour length. The curvature effects can be 
interpreted in terms of the polymer's backbone rigidity [71], etc. Upon Fourier- 
Laplace transform of the distribution function, Eq. (7.3.2), the propagator for 
the Klein-Gordon quantum field can be obtained [7, 72]. Hence, it is possible 
to provide an interpretation of the Ricci flow processes in terms of some dy- 
namics of quantum fields, etc. Whatever processes we may discuss, it should be 
clear that a specific physical process can be used only as some specific realiza- 
tion of the Ricci flow. Since such situation is common for all real life processes 
described by the equations of mathematical physics, the task lies in finding 
different physical situations in which this flow can be realized. For instance, 
since to conduct experiments on black holes [65] is unrealistic, one can think of 
analogous processes in liquid helium as described in the monograph by Volovik, 
Ref [73] . A great deal of additional condensed matter analogs of gravity-related 
phenomena can be found in Refs.[74] and [75] and references therein. 



8.2 Interplay between topology, geometry and physics 

To make our discussion complete, let us return back to the processes involving 
critical dynamics studied in this paper. In view of Theorem 7.20., the obtained 
scalar curvatures either go to zero (for expanding solitons) or blow up (for 
shrinking solitons) when time t — > oo [2]. Physically, the case of curvatures 
blowing up is not acceptable, however, e.g. see Eq.(3.1) and comments related 
to this equation 35 . Hence, physically acceptable manifolds must either be flat 
or hyperbolic in accord with Eq.(6.30) and Remark 6.8. Although this fact 
was discussed in our earlier works, [47,55], we would like to provide a sketch of 
relevant arguments in this subsection. 

We begin with two dimensional CFT models at criticality following Rcf[21], 
pages 340-344, and Ref. [47], Section 5, where additional details can be found. 
Many of these models can be obtained by some straightforward modifications of 
the simplest Gaussian model defined on the torus T 2 . If Wi and uj 2 are two toral 
periods, one can define their ratio r as r = — . The number r is necessarily 



complex with Imr = 



2 



= y . In terms of these notations, the dimcnsionlcss 



free energy T for the Gaussian model is obtained as follows: 

T= l -\ny\ V (t)\ 2 , (8.1) 

where the function r)(r) is the Dedekind eta function. Such a function is known 
as the modular function for the once punctured torus [76]. Thus, even though 



35 It may be of some relevance to black holes, however. This topic requires further study 
within the context of liquid helium and other type of condensed matter experiments just 
mentioned. 
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originally the torus topology for the Gaussian model was used, the actual com- 
putations for this model involving the use of the first Kronecker limit formula 
[47] leads to the final manifestly modular invariant result for T . The price for 
modular invariance is the switch in topology: from that for a flat torus (before 
the Kronecker limit is taken) to that for a punctured torus (after the limit is 
taken). The punctured torus can be obtained from some parallelogram of pe- 
riods in the complex plane C whose vertices are removed and whose sides are 
identified pair-wise. Since Euler's characteristic \ for torus is while that for 
the punctured torus is — l,the punctured torus represents an example of the 
hyperbolic surface (actually an orbifold [77]) as discussed in great detail in our 
work, Ref.[78]. Such a surface is not compact, however, while the results of the 
previous section were developed for compact surfaces. 36 In two dimensions the 
situation can be repaired by using the Schottky double construction widely used 
in designing of CFT models [79]. Use of this construction causes us to replace 
the punctured torus by the double torus, i.e. by the Riemann surface of genus 
2 which is also hyperbolic. 

As we had demonstrated in [47] (see also Remark 6.6.) even if one begins 
with the standard Euclidean cube (or any parallelepiped for that matter whose 
sides are pair-wise identified thus making it a T 3 ) use of the analog of the first 
Kronecker limit formula in 3d leads to a result similar to Eq.(8.1) indicating 
that in complete accordance with 2d case, the limiting manifold/orbifold is 
hyperbolic. This fact is consistent with results of our earlier work Ref. [55] on 
AdS-CFT correspondence where totally independent arguments were used to 
arrive at such a conclusion. 

We conclude this work with a brief discussion of how one can actually design 
compact hyperbolic 3 manifolds. A typical 3d hyperbolic manifold in which a 
3d CFT lives and evolves can be built using a 2d punctured torus. To do so, we 
would like to use the results of our earlier work, Ref. [78], in which we noticed 
that such a punctured torus is the Seifert surface for the figure eight (hyperbolic) 
knot. The evolution of surface automorphisms in fictitious time create a 3 man- 
ifold. A specific hyperbolic 3 manifold known as a 3-manifold fibered over the 
circle ( the puncture on the torus can be opened up so that it is homeomorphic 
to a circle S 1 ) as discussed in our work, Ref. [78]. The 3 manifold created in such 
a way will be cusped [55] and, hence, noncompact. To make such a manifold 
compact one can also use a hyperbolic double (analogous to a Schottky double 
as described in detail in Ref. [80] . Obviously, the same is achieved if one starts 
with a Schottky double torus and fibers it over the circle [80]. The direct 3d 
analog of the two dimensional punctured torus can be constructed as a product 
of the trice punctured sphere and an open interval (0,1) [80]. The obtained 
manifold is homeomorphic to a cube whose upper and lower faces along with all 
edges are removed. Such a designed 3 manifold is hyperbolic but unfortunately 
non compact. Clearly, Pcrelman's work is not limited to compact manifolds 
and, hence, his results still can be used. The only trouble with noncompact- 

36 This restriction is not essential however. It is possible to extend the theory of Ricci and 
Yamabc flows to the non compact surfaces [1,2]. Such an extension requires use of more 
sophisticated mathematical methods than those used in our paper. 
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ness lies in the fact that the Yamabe flow considered in our work leads to the 
Einstein-type spaces so that compactness is synonymous with being of Einstein 
type. We mentioned this already Ref.s [26,27] in which Yamabe functional for 
manifolds with boundary were discussed. More relevant to physical applications 
is the work by Mazzeo et al [81] and, also Ref.[82], in which the Yamabe prob- 
lem was studied for noncompact manifolds. It would be very interesting and 
challenging to study the Yamabe flow for noncompact manifolds discussed in 
Refs[81,82]. Many details of the construction of 3-manifolds and orbifolds are 
discussed in the exceptionally well written Ref. [77] . We conclude our paper by 
urging our readers to read this reference, which is also helpful for developing a 
solid understanding of the Poincare' and geometrization conjectures. 

Note added in proof. After this work was completed several impor- 
tant papers supporting and clarifying the results of our Section 8.2. came 
to our attention. In particular, in the paper by Long and Reid (Algebraic 
and Geometric Topology 2 (2002) 285-296) it is shown that all flat mani- 
folds can be looked upon as cusps of hyperbolic orbifolds, e.g. see our works, 
Ref.s [55], [78] for a quick introduction into cusps and Ref. [77] for introduction 
to hyperbolic orbifolds. Many additional details elaborating on the work by 
Long and Reid can be found in the PhD thesis by David Ben McReinolds 
(arxiv:math.GT/0606571). Some related material can also be found in the paper 
by X.Dai (arxiv:math.DG/0106172). 



42 



References 

[I] G.Pcrclman, The entropy formula for the Ricci flow and its geometric 
applications, arxiv:math.DG/0211159; ibid Ricci flow with surgery on 
three-manifolds, arxiv:math.DG/0303109; ibid Finite extinction time 
for the solutions to the Ricci flow on certain three- manifolds, 
arxiv:math.DG/0307245. 

[2] H-D.Cao, X-P.Zhu, A complete proof of the Poincare' and 

geometrization conjectures- Application of the Hamilton-Perclman 
theory of the Ricci flow, Asian Journal of Math.10 (2006) 165-492. 
[3] J.Morgan, G.Tian, Ricci flow and the Poincare' conjecture, 

arxiv:math.DG /0607607. 
[4] D.Fricdan, Nonlinear models in 2+e dimensions, 

Ann.Phys.163 (1985) 318-419. 
[5] R.Hamilton, Three-manifolds with positive Ricci curvature, 

J.Diff.Gcom. 17 (1982) 255-306. 
[6] B.Chow, P.Lu, Lei Ni, Hamilton's Ricci Flow, 

AMS Publishers, Providence, RI, 2006. 
[7] A.Polyakov, Gauge Fields and Strings, 

Harwood Academic Publ, New York, NY, 1987. 
[8] M.Green, J.Schwarz, E.Witten, Superstring Theory, Vol.1., 

Cambridge U. Press, Cambridge, UK, 1987. 
[9] D. Gross, J.Harvey, E.Martinec, R.Rohm, Heterotic string theory, 

Nucl.Phys.B267 (1986) 75-124. 
[10] B.Kleiner, J.Lott, Notes on Perelman's papers, arxiv: math.DG./ 
06056667. 

[II] L.Bessieres, Conjecture de Poincare' : la preuve dc R.Hamilton et 
G.Perelman, SMF-Gazette 106 (2005) 7-35. 

[12] M.Kapovich, Geometrization conjecture and the Ricci flow, 

www.math.ucdavis.edu/~kapovich/EPR/ricci.ps 
[13] B.Hatfield, Quantum Field Theory of Point Particles ans Strings, 

Addison- Wesley Publ.Co, New York, NY, 1992. 
[14] E.Abdalla, M.Abdalla, D.Dalmazi, A.Zadra, 2D-Gravity 

in Non-Critical Strings, Springer- Verlag, Berrlin,1994. 
[15] Yu.Nakayama, Liouvillc field theory- A decade after revolution, 

arxiv:hcp-th/0402009 
[16] R.Hamilton, The Ricci flow on surfaces, 

Contemp.Math.71 (1988) 237-262. 
[17] D.Amit, Field Theory, the Renormalization Group 

and Critical Phenomena, McGraw-Hill Inc., London, 1978. 
[18] C.Itzykson, J-B.Zuber, Quantum Field Theory, 

McGraw-Hill, London, 1980. 
[19] J.Lcc, T.Parkcr, The Yamabe problem, BAMS 17 (1987), 37-91. 
[20] K.Richardson, Critical points of the determinant of 

the Laplace operator, J.Funct. Anal. 122 (1994), 52-83. 
[21] P.Di Francesco, P.Mathieu, D.Senechal, Conformal Field Theory, 

Springer, Berlin, 1997. 



43 



[22] R.Wald, General Relativity, The University 

of Chicago Press, Chicago, 1984. 
[23] A.Kholodenko, K. Freed, Theta point ( "tricritical" ) region behavior 

for a polymer chain: Transition to collapse, J.Chem.Phys. 

80 (1984)900-924. 
[24] L.D. Landau, E.M.Lifshitz, Statistical Physics. Part 1. Course in 

Theoretical Physics. Vol.5., Pergamon Press, Oxford, 1982. 
[25] H. Bray, A. Neves, Classification of prime 3-manifolds 

with (7— invariant greater than RP 3 , 

Ann.Math. 159 (2004), 407-424. 
[26] J.Escobar, Conformal deformation of a Ricmannian metric 

to a scalar flat metric with constant mean curvature 

on the boundary Ann.Math. 136 (1992), 1-50. 
[27] K.Akutagawa, B. Botvinnik, Yamabe metrics on cylindrical 

manifolds, GAFA 13 (2003), 259-333. 
[28] R.M. Schoen, Variational theory for the total scalar 

curvature functional for Ricmannian metrics and related 

topics, LNM 1365 (1989), 120-184. 
[29] T. Aubin, Some Nonlinear Problems in Ricmannian Geometry, 

Springer, Berlin, 1998. 
[30] P.Dirac, General Theory of Relativity, 

Princeton U. Press, Princeton, 1996. 
[31] H. Yamabe, On a deformation of a Ricmannian structures on 

compact manifolds, Osaka Math. J. 12 (1960), 21-37. 
[32] A.Kholodenko, Quantum gravity and QCD in the light of works by 

Hamilton and Pcrclman, in preparation. 
[33] J.Schouten, Ricci-Calculus, Springer- Verlag, Berlin, 1954. 
[34] K.Belov, About extrinsic geometry of hypersurfaces, 

Soviet.Math.Izvestiya 104 (1971) 1-8 (in Russian). 
[35] C.Fedischenko, V.Chcrnyshcnko, About some generalization of 

spaces of constant curvature, Tensor Analysis and Appl. 11 

(1961) 269-276 (in Russian). 
[36] J. Zinn- Justin, Quantum Field Theory and Critical Phenomena, 

Clarendon Press, Oxford, 1989. 
[37] A.Polyakov, Quantum geometry of bosonic strings, 

Phys.Lctt.B 103 (1981), 207-210. 
[38] O.Alvarez, Theory of strings with boundary 

Nucl.Phys.B 216 (1983), 125-184. 
[39] W.Weisberger, Conformal invariants for determinants 
of Laplacians on Riemann surfaces, 
Comm.Math.Phys. 112 (1987), 633-638. 
[40] W.Weisberger, Normalization of the path integral measure 
and coupling constants for bosonic strings, 
Nucl.Phys. B 284 (1987), 171-200. 
[41] B.Osgood, R. Phillips, P. Sarnak, Extremals of determinants 

of Laplacians, J.Funct.Anal. 80 (1988), 148-211. 



44 



[42] B.Osgood, R. Phillips, P. Sarnak, Compact isospectral 

sets of surfaces, J.Funct. Analysis 80 (1988), 212-234. 
[43] B. Osgood, R. Phillips, P.Sarnak, Moduli space, heights and 

isospectral sets of plane domains, Ann. Math. 129 (1989), 293-362. 
[44] E. Onofri, On the positivity of the effective action in a theory 

of random surfaces, Comm.Math.Phys. 86 (1982), 321-326. 
[45] W. Beckner, Sharp Sobolev inequalities on the sphere and the 

Moser-Trudinger inequality, Ann.Math. 138 (1993), 213-242. 
[46] M.Goulian, M. Li, Correlation functions in Liouville theory. 

PRL 66 (1991), 2051-2055. 
[47] A.Kholodcnko, E.Ballard, From Ginzburg-Landau to Hilbert- 

Einstein via Yamabe, Physica A 380 (2007) 115-162. 
[48] S.Rosenberg, The Laplacian on a Riemannian Manifold. 

Cambridge Univ. Press, Cambridge, UK, 1997. 
[49] A.Besse, Einstein Manifolds, Springer- Verlag, Berlin, 1987. 
[50] T.Parker, S. Rosenberg, Invariants of conformal Laplacians, 

J.Diff.Geom. 25 (1987), 535-557. 
[51] K.Okikiolu, Critical metrics for the determinant 

of the Laplacian in odd dimensions, 

Ann.Math. 153 (2001), 471-531. 
[52] P.Chui, Height of flat tori, 

AMS Proceedings 125 (1997), 723-730. 
[53] S. Hawking, Zeta function regularization of path integrals in 

curved space, Comm.Math.Phys. 55 (1977), 133-148. 
[54] Y. Muto, On Einstein's metrics, 
J.Diff. Geom. 9 (1974), 521-530. 
[55] A.Kholodcnko, Boundary CFT, limit sets of Kleinian 

groups and holography, J. Geom. Physics 35 (2000) 193-238. 
[56] N.Koiso, Einstein metrics and complex structures, 

Inv.Math. 73 (1983), 71-106. 
[57] L. Landau, I.Khalatnikov, About anomalous sound 

attenuation near critical points of the phase transitions of 

second kind, Sov.Physics Doklady 96 (1954), 469-472. 
[58] P.Chaikin, T.Lubcnsky, Principles of Condensed Matter Physics, 

Cambridge U. Press, Cambridge, UK, 2000. 
[59] H.Schwctlick, M. Struwe, Convergence of the Yamabe flow 

for "large" energies, J. Reine und Angew.Math. 562 (2003), 59-100. 
[60] A.Isihara, Statistical Physics, Academic Press, New York, NY, 1971. 
[61] R.Miiller, Differential Harnac Inequalities and the Ricci Flow, 

European Math.Soc.Publ.House, Zurich, SW, 2006. 
[62] P.Topping, Lectures on the Ricci Flow, 

Cambridge U. Press, Cambridge, UK, 2006. 
[63] A.Borisenko, An introduction to Hamilton's and Pcrclman's work 
on conjectures of Poincare' and Thurston, 

www-mechmath.univer.kharkov.ua/geometry/preprint/perelman7.pdf 
[64] Y.Fujii, K-I.Maeda, The Scalar- Tensor Theory of Gravitation, 



45 



Cambridge U. Press, Cambridge, UK, 2003. 
[65] M.Headrick, T.Wiseman, Ricci flow and black holes, 

arxiv: hcp-th/0606086 
[66] H-D.Cao, Geometry of Ricci solitons, 

Chin.Ann.Math.27B(2) (2006) 121-142. 
[67] T.Ivey, Ricci solitons on compact three-manifolds, 

Diff.Geomctry and its Applications, 3 (1993) 301-307. 
[68] R.Hamilton, Non-singular solutions of the Ricci flow on 

three-manifolds, Comm. in Analysis and Geometry, 

Vol.7 (1999) 695-729. 
[69] P.De-Gennes, Scaling Concepts in Polymer Physics, 

Cornell Univ.Press, Itaca, NY, 1979. 
[70] R.Hamilton, Formation of singularities in the Ricci flow, 

Surveys in Diff.Geom.2 (1995) 7-136. 
[71] A.Kholodenko, Fermi-Bose transmutation: from 

semiflexible polymers to superstrings, 

Ann.Phys.202 (1990) 186-225. 
[72] A.Kholodenko, Statistical mechanics of deformable droplets 

on flat surfaces, J.Math.Phys.37 (1996) 1287-1313. 
[73] G.Volovik, The Universe in a Helium Droplet, 

Clarendon Press, Oxford, 2003. 
[74] W.Zurek, Cosmological experiments in condensed matter systems, 

Phys.Rcports 276 (1996) 177-221. 
[75] C.Barcelo, S.Libcrati, M.Visser, Analogue gravity from field 

theory normal modes, Class. Quantum Grav.18 (2001) 3595-3610. 
[76] H.Iwanicc, Topics in Classical Authomorphic Forms, 

AMS Publishers, Providence, RI, 1997. 
[77] D. Cooper, C. Hodgson, S.Kerhoff, Thrcc-Dimcnsional Orbifolds 

and Cone Manifolds, Japan Publ. Trading Co., Tokyo, 2000. 
[78] A.Kholodenko, Statistical mechanics of 2+1 gravity from 

Riemann Zeta function and Alexander polynomial: exact results, 

J.Gcom.Phys.38 (2001) 81-139. 
[79] C.Schweigerth, J.Fuchs, J.Walcher, Conformal field theory, 

boundary conditions and applications to string theory, 

arxiv: hep-th/0011109. 
[80] K.Matsuzaki, M.Taniguchi, Hyperbolic Manifolds and 

Klcinian Groups, Clarendon Press, Oxford, 1998. 
[81] R.Mazzeo, D. Pollack, K.Uhlenbeck, Moduli spaces of singular 

Yamabe metrics, AMS Journal 9 (1996) 303-344. 
[82] A.Chouikha, Ricci curvature and singularities of constant scalar 

curvature metrics, arxiv:math.DG/0501505. 



46 



